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PREFACE 


The  present  work  comprises  three  previously  published  papers  and  three  papers  which 
have  been  submitted  and  accepted  for  publication. 
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J M Young 
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W H Hook 

B R Bean 

(5)  Gjessing,  D T 
A G Kjelaas 

J Nordo 

(6)  Gjessing,  D T 
A G Kjelaas 
E Golton 


Tropospheric  structure  deduced  from  elec- 
tromagnetic and  acoustic  wave  propagation 
experiments.  Submitted  to  Geophysica  Nor- 
wegica. 

Study  of  divergence  in  the  boundary-layer 
using  optical  propagation  techniques, 
J Appl  Meteorol  13,  2,  242-8  (1974) 

Dispersion  and  spectra  of  gravity  waves 
probably  generated  by  a convective  storm, 
Tellus  27,  1 (1975).  In  press 


— Observation  of  waves  in  the  boundary-layer 
using  an  array  of  acoustic  sounders,  J At- 
mos Sci.  In  press 


— Spectral  measurements  and  atmospheric  sta- 
bility, J Atmos  Sci  26,  3,  462-7  (1969) 


Small-scale  atmospheric  structure  deduced 
from  measurements  of  temperature,  humi- 
dity and  refractive  index,  Boundary-Layer 
Meteorol  4,  475-92  (1973) 


' 

Paper  1 summarizes  the  potentials  of  different  electromagnetic  and  acoustic  wave 
propagation  experiments  as  diagnostic  tools  in  the  study  of  the  refractive-index,  tem- 
perature, humidity  and  velocity  fields  in  the  troposphere. 

Papers  2 and  3 discuss  measurements  of  horizontal  wind  velocity  using  optical  line- 
of-sight  propagation  techniques.  Paper  2 deals  with  measurements  of  horizontal  diver- 
gence while  in  Paper  3 the  spectrum  of  the  spatial  wind  together  with  the  spectrum 
of  air-pressure  are  used  to  deduce  the  propagation  speed  and  direction  of  progressive 
gravity  waves  in  the  troposphere. 

Paper  4 deals  with  the  potentials  of  an  array  of  vertically  pointing  acoustic  echo 
sounders  to  study  time- dependent  structures,  especially  waves,  in  the  boundary  layer. 

In  Paper  5 a forward  scatter  technique  (beam-swinging)  is  used  to  deduce  the  vertical 
spectra  of  refractive  index  and  air  velocity.  The  influence  of  atmospheric  stability  on 
these  spectra  is  discussed. 

ft  ' ~ 


I-  □ D, 


Paper  6 discusses  the  spatial  structure  of  refractive  index,  temperature  and  humidity 
in  the  lowest  2000  m of  the  troposphere  as  deduced  from  tethered  balloon  measure- 
ments. ^ 

With  the  exception  of  the  work  described  in  Papers  2,  3 and  4,  the  work  has  been 
carried  out  at  the  Norwegian  Defence  Research  Establishment. 

The  work  which  is  the  basis  for  Papers  2,  3 and  4 was  carried  out  during  a stay  as  a 
visiting  scientist  at  the  Wave  Propagation  Laboratory,  National  Oceanic  and  Atmos- 
pheric Administration,  Boulder,  Colorado,  and  was  sponsored  by  the  Royal  Nor- 
wegian Council  for  Scientific  and  Industrial  Research. 
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GENERAL  INTRODUCTION 

The  density  both  in  time  and  space  of  our  present  meteorological  observations  arc 
fairly  well  matched  to  the  synoptic  scale  weather  pattern. 

For  the  study  of  micro-  and  meso-scale  phenomena  it  is  evident  that  we  cannot  rely 
upon  ground-based  in  situ  measurements  and  twice  daily  radiosondes. 

As  it  is  neither  technologically  nor  economically  feasible  to  achieve  the  desired  reso- 
lution in  time  and  space  by  conventional  observation  techniques,  further  advances  in 
our  knowledge  about  the  atmosphere  at  micro-  and  meso-seales  are  strongly  depen- 
dent on  the  development  and  application  of  new  measuring  techniques. 

This  requirement  for  better  and  more  reliable  measurements  has  led  to  great  activity 
in  the  application  of  different  remote  sensing  techniques  which  are  based  on  the 
interaction  between  electromagnetic  or  acoustic  waves  and  refractive  index  fluctu- 
ations in  the  troposphere. 

The  progress  made  by  radio  communication  scientists  in  understanding  the  influence 
of  atmospheric  refractive  index  inhomogeneities  on  the  scattering  of  radiowaves,  has 
shown  that  radio  tropospheric  scatter  techniques  can  provide  important  information 
about  the  fine-scale  structure  of  the  troposphere. 

Highly  sensitive  radars,  mainly  developed  for  military  purposes,  have  proved  to  be 
extremely  valuable  tools  in  the  study  of  both  atmospheric  stratification  and  convec- 
tion. 

Based  on  the  experience  from  line-of-sight  communication,  a variety  of  optical  and 
microwave  methods  have  been  suggested  for  measuring  winds  and  refractive  index 
structures  in  the  boundary  layer. 

More  recently  it  has  also  been  demonstrated  that  acoustic  echo  sounding  can  be  a 
valuable  tool  both  in  the  study  of  wind  in  the  troposphere  and  for  monitoring 
temperature  inversion  structures  and  convection. 

The  present  work  describes  different  applications  of  linc-of-sight  optical  and  micro 
wave  propagation,  Doppler  radars,  forward  scatter,  and  acoustic  echo  sounding  in  the 
study  of  the  tropospheric  refractive  index,  temperature,  humidity  add  velocity  fields. 

It  is  shown  that  carefully  designed  experiments  involving  different  kinds  of  electro- 
magnetic and  acoustic  wave  propagation  can  give  reliable  and  very  useful  information 
regarding  the  structure  of  the  troposphere. 
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TROPOSPHERIC  STRUCTURE  DEDUCED  FROM  ELECTROMAGNETIC  AND 
ACOUSTIC  WAVE  PROPAGATION  EXPERIMENTS 


SUMMARY 

I his  paper  reviews  briefly  the  potentials  of  the  atmospheric  remote  sensing  techniques 
where  the  refractive  index  inhomogenities  play  the  dominant  role  in  producing  fluctu- 
ations in  the  received  signal.  The  techniques  discussed  are  optical  and  microwave 
line-of-sight  propagation,  acoustic  echo  sounding,  and  microwave  back-  and  forward 
scattering. 

Methods  for  measurements  of  variability  of  wind,  temperature  and  water  vapour  arc 
discussed  and  emphasis  is  laid  on  methods  which  have  been  proven  through  field 
experiments. 

lhe  theories  of  line-of-sight  propagation  and  scattering  are  also  reviewed  briefly. 


INTRODUCTION 

Further  advances  in  our  knowledge  about  the  atmosphere  are  strongly  dependent  on  the 
development  and  application  of  new  measuring  techniques.  The  shortcomings  of  point 
measurements  and  radiosondes  have  become  more  and  more  evident. 

For  this  reason  there  has  been  a rapid  development  of  a series  ot  new  techniques  during 
the  last  decade.  The  common  feature  of  one  group  of  these  techniques  is  that  they 
transmit  an  electromagnetic  or  acoustic  wave,  which  interacts  with  the  atmosphere.  A 
portion  of  this  energy  is  scattered  and  detected  at  the  receiver.  Information  about  the 
atmosphere  is  then  derived  from  the  characteristics  of  the  received  signal.  The  advan- 
tages of  such  a measuring  technique  are  many  (Little,  1973): 

— The  measuring  system  does  not  significantly  modify  (or  interfere  with)  the  medium 
being  measured 

— The  measurements  are  taken  remotely,  i e no  (in  situ)  instrumentation  has  to  be 
brought  to  the  region  of  the  atmosphere  to  be  measured 

— The  relevant  atmospheric  parameters  are  measured  in  one,  two  or  three  dimensions, 
mostly  with  excellent  resolution  in  both  time  and  space 

Further,  this  type  of  measurement  may  provide  the  following  information  about  the 
atmospheric  structure: 

a)  line  integral  the  integrated  value  of  a parameter  along  a line  through  the 

atmosphere 

b)  line  average  the  average  value  of  the  parameter  along  some  line  of 

known  length 

c)  line  profile  the  distribution  of  the  parameter  along  some  line  of  known 

length 

d)  2-dimensional  coverage  the  distribution  of  the  parameter  over  a plane 

e)  3-dimensional  coverage  the  distribution  of  the  parameter  in  space 

f)  spectrum  the  spatial  power  spectrum  of  a parameter 

g)  flux  the  rate  at  which  mass,  momentum  or  heat  is  transported 
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The  success  of  the  measuring  technique  is  dependent  on  the  strength  and  type  of  the 
interaction  of  the  wave  with  the  atmosphere.  We  may  distinguish  between: 

a)  a spectroscopic  interaction,  in  which  the  radiation  reacts  with  the  molecules,  i e an 
absorption  and  reradiation 

b)  a refractive  index  interaction,  where  the  scattering  of  the  incident  wave  is  caused 
by  inhomogenities  in  the  refractive  index  field,  or,  more  specifically,  the  three- 
dimensional  spectrum  of  refractive  index  fluctuations 

For  this  short  report  we  shall  only  consider  measuring  techniques  in  which  refractive 
index  inhomogenities  play  the  dominant  role.  The  discussion  is  also  narrowed  to 
’’sensors”  which  are  located  on  or  near  the  ground. 

The  techniques  discussed  are  line-of-sight  optical  and  microwave  propagation,  acoustic 
echo  sounding,  pulse  Doppler  and  FM— CW  radars, and  forward  scatter. 

Scattering  usually  implies  that  no  signal  will  be  received  in  the  absence  of  refractive 
index  inhomogenities,  while  in  the  case  of  line-of-sight  propagation  the  field  at  any 
point  will  be  the  sum  of  the  direct  wave  and  the  many  scattered  waves,  scattered  and 
refracted  by  local  variations  in  the  refractive  index  field. 

Each  technique  is  discussed  separately  and  emphasis  is  laid  on  capabilities  which  have 
been  demonstrated  through  field  experiments. 


2 THE  THEORETICAL  BACKGROUND 

In  this  section  we  shall  derive  the  necessary  equations  governing  EM  and  acoustic  wave 
propagation  in  the  atmosphere. 

2.1  Refractivity  and  refractive  index 

The  radio  refractivity  N for  air  is  given  in  terms  of  refractive  index  n as 

N = (n-l)-lO2 * * * 6  = ^ (p+4810  ^)  (2.1) 

where 

T — temperature  in  degrees  Kelvin 
p — total  pressure  in  mb 
e — partial  water  vapour  pressure  in  mb 

For  acoustic  waves  the  refractivity  Na  is  given  by 

Na  = (na— 1)-106  * 106  _ 0.14  J ] (2.2) 

where 


cQ=  20.5  T01/2  is  a reference  sound  velocity  at  a reference  temperature  TQ 
0 — angle  between  the  wave  and  wind  vector 
V — wind  speed 
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1 hr  optical  refractive  index  nu  can  l>c  written  as 
Nu  - (no-l)10*  ^ 79  2 


(2.3) 


MASNITUOi 
CHANS IS 

CHANGES  IN  REFRACTIVITY  IN  NUNITS 

RADIO 

OPTICAL 

ACOUSTIC 

AT  • 1 It 

1 

1 

1.7 -ID’ 

AV  • 1 KM  1 

- 

- 

Ml* 

A*  * 1 Mb 

4 

4 ir» 

1.4- ISJ 

Table  2.  I Helat ire  ralucs  of  Cn  for  acoustic  radio 

and  optical  wares 


The  influence  of  small  changes  in 
the  atmospheric  parameters  on  the 
different  refractive  indexes  is 
shown  in  Table  2.1.  Thus,  for 
example,  acoustic  waves  are  more 
than  three  orders  of  magnitude 
more  sensitive  to  changes  in  tem- 
perature than  either  radio  or  opti- 
cal waves.  While  acoustic  waves 
are  sensitive  to  wind  fluctuations, 
the  refractive  index  for  electro- 
magnetic, including  optical,  waves 
is  independent  of  wind. 


I he  refractive  itulex  can  lx-  decomposed  into  a mean  n plus  a fluctuating  component  n( 
as 


n(r)  = n(rj+n,(r)  (2.4) 

Often  it  is  very  useful  to  represent  the  fluctuating  component  in  terms  of  a covariance 
function  and  its  associated  spectrum  or  in  terms  of  a stochastic  Fourier  Stieltjes 
integral. 

The  covariance  function  is  defined  as 

C(r)  « n,(  r, ) n,(r,  + r)  (2.5) 

and  its  associated  spectrum  as 

♦nffc)  = — ///  C(r)exp(  iiv  r)dr  (2.6) 

In  terms  of  a stochastic  Fourier  Stieltjes  integral,  n,  (r)  can  be  represented  as 

ni  (r)  * f/J  A(lv)exp(iKr)dK  (2.7) 


and  the  spectrum  d>n(K)  can  be  obtained  by  evaluating  A(K).V  (K), where  the  asterisk 
denotes  the  complex  conjugate. 


2.2  Propagation  of  radio  and  optical  waves 

We  consider  a medium  where  the  refractive  index  is  expressed  by  equation  (2.4). 

I he  propagation  of  an  electromagnetic  wave  in  this  medium  is  governed  mathematically 
by  Maxwell’s  equations,  which  can  be  combined  to  give  the  vector  wave  equation 

V:£+kV(r)K  + 2V  £.V(n(r))  = 0 

where  it  is  assumed  that  the  electric  field  F.  is  sinusoidal,  and  k is  the  wavenumber. 


(2.8) 
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l he  last  term  in  hqu.it urn  (2.8)  describes  the  polarization  of  the  field.  If  the  scattered 
depolarized  field  is  much  smaller  than  the  field  in  the  original  direction  of  polarization, 
hquation  (2.8)  reduces  to  the  scalar  wave  equation 

7JE+k1n,(r)  E = 0 (2.9) 

Different  approaches  have  been  used  to  obtain  an  approximate  solution  of  this  equation 
when  the  medium  is  weakly  inhomogeneous: 

Geometrical  or  ray  optics 

Method  of  small  perturbations  or  the  first  Horn  approximation 
Method  of  smooth  or  slowly  varying  perturbations  (the  Kvto\  method) 

For  details,  the  reader  is  referred  to  Strohbehn  ( l‘H>8). 

In  the  case  of  large  perturbations,  ter  hniques  for  solving  the  equation  have  not  vet  been 
satisfactorily  developed. 


Scattering  of  elei  tromagnetic  waves 


Consider  a volume  element  d'r  within  a scatter  volume  V filled  with  a random  medium 
with  a refractive  index  as  defined  in  Equation  (2.4).  It  illuminated  hv  an  incident  field 
Eu  the  element  becomes  polarized,  giving  rise  to  a dipole  moment  dP-  n,  F.  d’r.  At  a 
distance  R the  dipole  moment  results  in  a polarization  potential.  If  R V1  ’ the 
scattered  field  Es  resulting  from  all  the  elementary  scattering  elements  can  be  expressed 
.is 


F.  ks  K°  A(K) 
s <tr  R 


(2.10) 


0 and  ks  are  the  ineident  and  scattered  wave  vei  tors.  K k()  ky 
A(K)  is  obtained  from  the  inverse  transform  of  Equation  (2.7). 
k = | k0l  I ksl  = 2>r/X,  where  X is  the  wavelength. 

R is  the  distance  from  the  scattering  volume  to  the  observation  point  at  which  the  field 
is  equal  to  F.s.  Figure  2. 1 shows  the  scattering  geometry  . 


Figure  / Scatter  geometry' 


Die  magnitude  of  the  scattering 
vector  K kt)  ks  is 

K 2k  sin  J (2.11) 

where  0 is  the  angle  between  k 
1 r*  o 

and  ks,  i e the  scattenng  angle. 

l'he  scattering  cross-section  a 
which  is  a measure  of  the  fraction 
of  the  incident  power  scattered  in 
a given  direction,  can  be  obtained 
\ M by  evaluating  F.s-F.s"  (G jessing, 
1962).  In  terms  of  the  refractive 
index  fluctuations  the  following 
expression  of  the  scattering  cross- 
section  can  be  given 


where  X — polarization  angle 

'^n^)  three-dimensional  power  spectrum  of  the  refractive  index  fluctu- 
ations tt) 

from  F.quation  (2.12)  we  observe  that  the  intensity  of  the  scattered  wave  for  a given 
geometry  and  a given  wavelength  is  determined  by  one  spectral  component  in  the 
three-dimensional  refractive  index  spectrum,  corresponding  to  a scale  t of 


which  is  Bragg’s  condition  for  diffraction  from  a spatial  structure. 

for  pure  backscatter  (0  = 1 80  ) the  intensity  of  the  scattering  is  determined  by  the 
refractive  index  inhomogenities  of  a scale  size  corresponding  to  X/2. 

2.2.2  Iine-of-sight  propagation,  radio  and  optical 

1 he  scalar  wave  equation  can  be  solved  by  one  of  the  techniques  listed  in  section  2.2. 
Instead  of  using  one  of  these  approaches,  we  will  apply  a method  described  by  Lee  and 
Harp  (1969).  Ihis  method,  which  is  based  on  phase  screens,  leads  to  an  easier  physical 
understanding  of  the  problem  than  do  the  more  mathematical  approaches. 

In  line-of-sight  propagation  it  is  useful  to  decompose  the  refractive  index  field  n,  into  a 
plane  normal  to  the  propagation  direction  z,  expressed  in  terms  of  a two-dimensional 
stochastic  Fourier— Stieltjes  integral,  i e 

"I  (x.y.z)  = //  exp(iK-/J)  t>(z,K)d2K  (2.14) 

where  p is  a position  vector  in  this  plane,  /?=(x,y)  and  K=(KX.K  ).  f(z,K)  is  the  Fourier 
component  in  the  same  plane. 

Ihis  Fourier  component  will  cause  a small  phase  change  in  the  incident  wave  as  it  passes 
through  the  plane.  Ihis  phase  perturbation  produces  perturbation  in  the  phase  of  the 
received  field,  lhe  field  in  the  receiver  plane  due  to  all  the  Fourier  components  of  the 
refractive  index  and  all  of  the  layers  of  thickness  dz  between  transmitter  and  receiver  is 
then  for  a sphtneal  wave  (Lee  and  Harp,  1969) 

E = 1— ik  / dz  // t;(z,K)  exp(i  ^|^=^]Cos(K  /5^)d2K  (2.15) 

where  L is  the  length  of  the  path. 

Assuming  E = A exp(iS)  where  A = 1 +A|  and  S = kL+S|  we  can  express  the  amplitude 
and  phase  fluctuations  in  the  receiving  plane  in  terms  of  spatial  covariance  functions  as 
defined  in  Equation  (2.5). 

The  spatial  covariance  functions  can  also  be  derived  by  evaluating  A,  A,"  and  S,  S,  * 
(Tatarski,  1961). 

Assuming  that  the  refractive  fluctuations  arc  statistically  homogeneous,  the  two-dimen- 
sional spectrum  F(z,K)  of  the  refractive  index  field  can  be  obtained  by  evaluating 
v(z,K)  ^(z,K). 
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For  the  case  of  isotropy,  the  two-dimensional  spectrum  is  related  to  the  covariance 
function  as 


h‘(*.K)  - /qz,p)J0(Kp)pdp  (2.16) 

where  J0  is  a Bessel  function  of  zero  order. 

.Vssuming  isotropy,  the  spatial  covariance  function  of  aiiipliuule  and  phase  can  then  be 
expressed  in  terms  of  the  three-dimensional  refractive  index  spectrum  as 

Cs(p)  = 4irkJ  / / K4>  (K)  cos|  I j0(^)  dKdz  (2.17) 

o o JkL  w L 

c:a(p)  = 4 wI  2 k 2 //  K4>  (K)sin|!^Uc-il|  J0(^fi?)dKdz  (2.18) 

o o JkL  I. 

If  the  refractive  index  fluctuations  are  homogeneous  over  the  path,  the  integration  over 
z may  Ik-  carried  out. 

I ruler  the  assumption  of  Iaylor  s hypothesis,  the  temporal  power  spectrum  of  the  phase 
and  amplitude  fluctuations  can  be  obtained  by  substituting  into  Equations  (2.17)  and 
(2.18)  p = vxr,  where  vx  is  the  cross-wind  component,  and  taking  the  Fourier  transform 
of  Cs(vxr)  and  Ca(vxr).  r is  the  time. 


Propagation  of  acoustic  waves 


Similar  to  the  discussion  of  propagation  of  radio  and  optical  waves,  we  consider  a 
turbulent  medium  where  the  refractive  index  is  composed  of  a mean,  n , plus  a 
fluctuating  component  na  a 


na  = na+na,  (2.19) 

where  the  fluctuating  part  na  may  be  written  as 
u I 1 T, 

n‘*t"~^~2T^  (2-20) 

° o 

when  I r»a , I « |ffal 

where  cQ  is  the  velocity  of  sound,  corresponding  to  the  average  temperature  T . T,  is 
the  temperature  fluctuation  from  this  average  value  and  u(  is  the  perturbation  velocity 
of  zero  mean. 


I he  propagation  of  sound  in  this  medium  is  described  by  the  equations  of  hvdrotiyna- 
mics,  from  which  a perturbated  Helmholtz  equation  analogous  to  Equation  (2.8) 
(Clifford  and  Brown,  1970;  Nlonin,  1962)  can  be  written  as 

V’Ps+k’V  2k’nai  P0-2ikP0  6^I  = 0 (2.21) 

where  Ps  — scattered  acoustic  field  (acoustic  pressure) 

PG  = A0c‘kz  — incident  field 

I he  fourth  term  in  Equation  (2.21)  is  analogous  to  the  polarization  term  in  EM 

propagation. 


II 


While  the  polarization  term  was  very  small  in  the  EM  case,  it  has  to  be  included  in 
acoustic  propagation  since  the  acoustic  wave  is  essentially  longitudinally  polarized  and 
accordingly  will  be  affected  by  wave  motion  parallel  to  gradients  of  refractivity. 


2.3. 1 Scattering  of  acoustic  waves 


I he  scattering  of  sound  waves  in  a turbulent  medium  is  analogous  to  scattering  of 
electromagnetic  waves.  The  scattered  acoustic  lield  Ps  can  be  expressed  in  terms  of  the 
three-dimensional  Fourier  transform  of  the  refractive  index  as 


Ps  = 


k2P 

*sro 

2rrr 


!na,(r)- 


fin 


ai 

~sr 


1 e'^'^r 


(2.22) 


I he  scattering  cross-section  is  found  analogously  to  the  EM  case  by  evaluating  Ps-Ps*.  In 
terms  of  the  refractive  index  fluctuations  the  following  expression  of  the  scattering 
cross-section  can  be  given  (Tatarski,  1971) 


7T 

2 


cos2  I *V(K) 


(2.23) 


with  the  same  definition  of  a as  in  the  EM  case. 


^(K)  and  4>y(K)  are  the  three-dimensional  power  spectral  densities  of  temperature  and 
wind  fluctuation,  respectively.  The  expression  for  K is  similar  to  the  EM  case,  i c,  K = 
2k  sin  0/2,  where  k is  the  wavenumber  of  the  sound,  and  T0  is  the  temperature 
corresponding  to  the  sound  velocity  c0>  From  Filiation  (2.23)  we  sec  that  the  back- 
scatter  (0=180°)  is  produced  by  temperature  inhomogenities  since  the  factor  cos2  0/2  is 
zero.  For  forward  scatter  we  have  contributions  from  both  temperature  fluctuations  and 
turbulence.  I hat  the  velocity  fluctuations  do  not  contribute  to  the  backscatter,  but  only 
in  the  forward  scatter  case,  is  a direct  consequence  of  the  incompressibility  of  the 
turbulent  flow. 


2.4  Doppler  techniques 

For  both  EM  and  acoustic  wave  propagation,  Doppler  techniques  can  be  used.  With  a 
scatter  geometry  as  shown  in  Figure  2.1,  the  Doppler  shift  is  given  by 

Af=2?(l^s-0V  (2.24) 

or  (Ottersten  et  al,  1973) 

Af  * ^ sin  ^ cos(3V 

if  we  neglect  the  difference  in  wavelength  between  the  incident  and  scattered  wave. 

Hie  angle  (3  is  the  angle  between  the  wind  vector  V and  dir  (ks—kQ).  In  the  case  of 
backscatter  (0/2  = 90),  the  Doppler  shift  is  a measure  of  the  radial  velocity. 

Since  the  obtained  velocity  is  not  observed  at  a point  in  space  but  over  a volume,  we  will 
observe  a Doppler  spread  or  a Doppler  spectrum  due  to  fluctuations  in  the  wind  velocity 
throughout  the  volume  and  also  from  the  fact  that  i^-jc  varies  through  the  volume  (for 
forward  scatter).  I he  first  moment  of  this  spectrum  gives  the  average  wind  velocity, 
while  the  variance  of  the  wind  velocity  can  be  derived  from  the  second  moment. 

My  scanning  a radar  in  azimuth,  several  radial  components  of  the  wind  can  be  obtained. 
This  technique,  which  is  called  VAI)  (Velocity— Azimuth  Display),  was  first  proposed  by 
Lhermittc  and  Atlas  (1961).  The  geometry  of  the  VAI)  technique  is  shown  in 
Figure  2.2. 


Figure  2.2  Geometry  oj  the  VAD  technique  (from  l. hermit te  and  Atlas.  1961) 


The  mean  wind  direction  and  speed  as  obtained  from  the  VA1)  technique  can  be  written 
as 


< 

X 

II 

V„  cos  (P-Pa)  + Vf  sin  0 

where 

Vr 

— radial  velocity 

VA 

— horizontal  wind 

6 

— elevation  of  the  scan 

Vf  sin  0 

--  DC  component  due  to  vertical  wind 

Browning  and  VVexler  ( 1968)  have  shown  that  by  expanding  the  radial  velocity  V^(a)  as 
a function  of  azimuth  in  Fourier  scries 


Vk(a)  = AM+ 


2 

E (Ak  cos  ka  + Bk  sin  ka) 


(2.26) 


the  Fourier  coefficients  are  related  to  the  mean  wind  components  u v , to  the 
horizontal  divergence  3u/3x  + 3v/3y,  to  the  stretching  deformation  3u/3x  dv/dy  and 
to  the  shearing  deformation  3u/3y  + 3v/3x  as 


v„  = A | /cos  0 


uD  = Hi  /cos  0 


3u  + 3v  _ 
3x  3y 


2 

Rcos  30 


( Aq  — Vpsin  0) 


(2-27) 


3u  _ 3v  _ _ 2A2 

3x  3y  Rcos20 

3u  3v~  - 2Bt 

3y  3x  Rcos  20 

If  we  want  to  determine  the  three-dimensional  wind,  three  radars  at  different  locations 
observing  the  same  volume  of  the  air  arc  nee  led. 


IJicrmittc  and  Miller  (1971)  have  proposed  a system  where  two  Doppler  radars  arc 
scanned  in  such  a way  that  the  two  radar  beams  stay  in  the  same  plane  during  a given 
time  intcrv;d.  I'his  method,  which  they  called  t'OPLAN  scanning,  gives  the  two-dimen- 
sional wind  in  a horizontal  plane. 


r 
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2.5  Description  of  ihr  atmosphere  in  terms  of  EM  and  acoustic  wave  propagation  para- 
meters 

1 he  influence  of  the  atmosphere  enters  into  the  equation  for  scatter-  and  line-of-sight 
propagation  of  electromagnetic  waves  through  the  refractive  index  spectrum  ‘J'^K)  and 
the  cross-path  wind  component  vx.  For  acoustic  waves,  the  scattering  is  determined  by 
the  temperature  spectrum  ‘I* ( K ) and  the  turbulence  spectrum  «J*\.(K). 

Assuming  isotropy  ami  incompressible  turbulence,  these  spectra  are  usually  expressed  as 
(Tatarski,  1961) 


. 


<1>,.(K)  = 0.033  C,  2 K "J 

(2.28a) 

0>n(K)  = 0.033  Cn2  K 1,3 

(2.28b) 

«I*V(K)  O.Otil  Cv2  K " 1 

(2.28c) 

where  Cj  2,  Cn2  and  Cy2  are  the  so-called  structure 
live  index  and  wind,  respectively. 

constants  for  temperature,  refrac- 

The  structure  constants  are  a measure  of  the  intensity  of  the  fluctuations  and  they  are 
defined  by 

C, 2 =1  T(x)-T(xfr)  |3 

1 r,,J 

(2.29a) 

c 2 = lyLll^x+rl  I2 

» 1 rl/J  1 

(2.29b) 

and 

r , 2 = I »jxl  u(xTr)  . 3 
' rl/3 

(2.29c) 

u is  the  wind  speed  and  r is  a distance  measured  along  the  x-axis.  The  x-axis  may  be 
chosen  in  any  arbitrary  direction  to  suit  the  experiment.  The  structure  constant  can  also 
he  expressed  in  terms  of  the  energy  dissipation  f as  (Tatarski,  1961) 

Cv2  = a2e2/3  (2.30) 

where  a is  a constant. Equations  (2.28a,h,c)  are  valid  within  some  range  of  K which  is 
called  the  inertial  subrange.  The  inertial  subrange  is  defined  by  an  outer  scale  L and  an 
inner  scale  CQ,  shown  schematically  in  Figure  2.3. 
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Figure  2.3  .1  schematic  diagram  of  the  inertial 
subrange 


The  inner  scale  L'u  which  is  where 
the  dissipation  starts,  is  in  the 
atmosphere  assumed  to  be  of 
order  mm,  while  the  outer  scale 
L0  which  is  the  scale  of  energy 
input,  is  for  isotropic  turbulence 
assumed  to  be  comparable  with 
the  height  above  ground,  but  very 
little  information  exists  on  the 
outer  scale. 


If  the  magnitude  of  the  wave- 
number vector  lies  inside  the  iner- 
tial subrange,  i e 


2jt 


< 2k  sin 


0 <-  2ff 
5 «o 


A 


then  the  expression  for  the  different  spectra  given  in  Equation  (2.28)  can  l>e  suhstitnted 
in  the  expressions  describing  the  scatter  and  line  of  sight  propagation. 


3 THE  POTENTIALS  OF  WAVE  PROPAGATION  EXPERIMENTS  IN  DERIVING 
INFORMATION  ABOUT  TIIE  ATMOSPHERE 

In  the  previous  set  lion  we  discussed  briefly  the  theorem  al  framework  for  both  EM  and 
ac  oust  it  wave  propagation.  We  observed  that  the  alinospherh  stun  lures  whnli  enter  the 
different  propagation  equations  are: 

the  refractive  index  spet  tnim  »l>  (k) 
the  temperature  spectrum  'I',  (k),  and 
the  turbulence  spectrum  'Ey  (k),  in  addition  to 
the  cross  beam  component  of  mean  wind  velocity 

In  this  section  we  shall  discuss  the  diagnostic  potentials  of  each  of  the  techniques 
described  in  the  previous  section.  • 

3.1  l.ine-of-sight  propagation 

We  have  shown  that  line  of  sight  propagation  is  alio  ted  b\  the  i loss  beam  wind  veloi  ilv 
and  by  the  refractive  index  spectrum.  These  far  too  will  now  be  disi  ussed  m some  detail. 

3.1.1  Determination  of  the  refrai  live  index  spectrum  ‘I*  (k) 

The  observed  fluctuations  in  amplitude  and  phase  of  a radio  wave  or  an  optical  wave 
over  a line  of  sight  path  are  determined  by  the  refrai  tne  index  spei  Hum 

The  intensity  of  the  refractive  index  fluctuations,  l^3,  averaged  over  the  transmission 
path  can  be  obtained  by  measuring  the  amplitude  of  the  received  signal.  This  line 
average  of  G 3 is  a weighted  average  weighted  both  along  the  path  and  in  the 
wavenumber  k.  The  amplitude  fluctuations  are  most  sensitive  to  stale  si/es  correspon 
ding  to  the  first  Fresnel  zone  XI-  Phase  fluctuations  are  sensitive  to  all  scale  sizes. 

But  since  most  of  the  energy  in  the  spectrum  is  associated  with  wavenumbers 
K * 2jt/I.0,  phase  fluctuations  are  mostly  determined  by  the  ontei  scale. 

A spatial  distribution  of  l',,3  along  the  path  may  be  taken  into  consideration  by 
assuming  a series  of  step  functions  along  the  path  (Harp,  I ‘>71).  Any  non  uniformity  in 
the  form  of  the  refractive  index  spectrum,  however,  will  make  tire  interpretation  verv 
difficult. 

The  form  of  <t»n(K)  averaged  over  the  path  can  Ire  obtained  by  solving  tire  spatial 
covariance  function  in  Equations  (2.17)  and  (2.18).  In  some  cases  tins  can  easily  be 
done,  but  in  most  cases  we  have  to  use  numerical  techniques  because  we  end  up  with 
non-linear  integral  equations.  For  details  the  reader  is  referred  to  Strolrbelrn  (1973)  and 
Lee  (1969). 

The  influence  on  line  of  sight  propagation  of  ‘h|((kl  Iras  been  examined  bv  Strohbehn 
(1970)  and  lx-e  and  Harp  ( 1 969).  Figure  3. 1,  which  is  taken  from  l.ee  and  Harp  ( 1 969), 
shows  the  influence  of  the  form  of  the  spectrum  on  the  spatial  amplitude  covariance 
function.  They  found  that  the  amplitude  covariance  function  was  reasonably  sensitive  t» 
the  exponent  n (spectrum  slope)  when  the  spectrum  was  written  as 


«*»n(K)  - k " 


Any  non  uniformity  along  the  path  ol  the  refractive  index  spectrum  will  again  complicate 
the  picture  and  we  may  conclude  th.it  linc-of  sight  propagation  is  not  specifically  well 
suited  for  determining  the  refractive  index  spectmm. 


hf^ire  3.  I I mphtude  covanam  c June turns  for  different  spec  fra  4*(  K ) with  experimental 
data  ( I 2S  km.  k 7 Hi  m 1 ) (from  l ee  and  llarp.  1969) 

Notr  that  u main  svavriiuiiitwr. 


Determination  of  wind  speed  from  line-of-sight  experiments 

( ontrary  to  the  case  of  determining  the  refractive  index  spectrum,  line-of-sight  seems  to 
Ik1  a very  useful  sensor  in  determining  the  average  cross  path  wind  speed. 

Two  different  methods  have  been  used: 

a)  I he  slope  method  discussed  hy  Shcn  ( I D 70)  and  l-awrcnce  et  al  ( 1 0 7 ) . 

I his  method  uses  the  slope  at  zero  lag  of  the  lagged  cross-covariance  function  of 
the  scintillation  at  two  points  in  a receiving  plane  normal  to  the  path  as  a measure 
of  the  wind  speed  across  the  path. 

h)  l he  delay  method  which  uses  the  time  lag  corresponding  to  the  maximum  of  the 
same  cross-covariance  function. 

1 he  concept  of  these  techniques  is  schematically  shown  in  Figure  3. '2-  Figure  .S..H  shows 
a comparison  between  the  transverse  wind  obtained  from  the  slope  method  and  a simple 
average  from  seven  propellor  type  anemometers  spaced  along  the  path.  By  using  three 
such  systems  arranged  to  form  a triangle,  horizontal  divergence  has  been  successfully 
measured  (Kjelaas  and  Ochs,  D>74).  The  same  technique  has  been  used  in  a study  of  the 
spatial  horizontal  wind  in  connection  with  gravity  wave  induced  wind  fluctuations 
(Kjelaas  et  al,  1 074b). 


Figure  3.2 

.1  schematic  diagram  of  the  normal 
ized  covariance  function,  showing 
delay  to  the  peak  and  the  slope  at 
zero  lag.  While  the  former  is  fre 
quently  used  to  measure  pattern 
drift  velocities,  the  latter  is  more 
suitable  in  the  optii  al  case  where  the 
pattern  decays  rapidly  (from  law 
rence  et  al,  1972) 
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hgurt'  .>.  i .1  comparison  of  the  average  witui  speed  ( 1',,)  measured  optically,  with  the 
average  of  the  readings  of  six  CM  anemometers  ( l'(1)  uniformly  spaced  along 
the  optical  path  (from  Kjelaas  and  Ochs,  197-4) 

Both  the  anemometer*  anil  the  optical  system  measure  the  component  of  wind  which  is 
horizontal  and  at  right  angles  to  the  path.  The  path  is  S00  m long,  oriented  north  south, 
and  S m above  the  ground. 


Measurement  of  vertical  wind  profiles  using  scintillations  of  light  from  stars  has  been 
attempted,  but  this  technique  is  not  sufficiently  well  developed  to  give  detailed  results. 
It  probably  requires  detailed  information  about  the  refractive  index  spectrum  as  a 
function  of  height,  which  is  difficult  to  obtain,  Ciossard  ( 19(10)  and  Gjessing  et  al  ( 1 97.’!) 
have  found  that  the  form  of  the  spectrum  varies  in  a complicated  manner  with  height. 

S.2  Radars 

According  to  Kquation  (2.12)  the  backseat tcred  signal  is  determined  by  the  refractive 
index  spectrum,  i e by  the  spectral  density  at  the  wave  number  corresponding  to  one 
half  wavelength,  K = 4tr/\.  The  direction  of  K is  in  the  radar  radial  direction. 

In  the  study  of  the  backscattcr  front  refractive  index  fluctuations,  two  different  regimes 
of  atmospheric  turbulence  are  distinguished:  the  convection  regime  and  the  stable 
regime. 

3.2.1  Detection  of  clear  air  convection 


i 


In  the  case  of  convection  the  radar  backseat  ter  is  produced  by  refractive  index  fluctua- 
tions associated  with  rising,  warm,  moist  air  produced  by  surface  heating.  Positive 
buoyancy  forces  contribute  to  the  generation  of  small-scale  turbulence. 


r> 


If  the  wavenumber  K = 4rr/X  falls  within  the  inertial  subrange,  the  radar  cross-section  for 
a region  of  refractive  index  fluctuations  will  be  proportional  to  Cn2  averaged  over  a 
volume  determined  by  the  beamwidth  and  the  pulse  length.  As  the  structure  constant 
Cn2  is  connected  with  dissipation  of  energy  and  refractivity  fluctuations  through 
Equation  (2.30),  radar  measurements  may  provide  a complete  description  of  the 
microstructure  in  the  inertial  subrange.  For  more  details,  the  reader  is  referred  to 
Ottersten  (1969). 

Atlas  and  Hardy  (1966)  and  Konrad  (1970),  among  many  others,  have  observed  free 
convection  using  highly  sensitive  radars. 

Recently  Noonkester  et  al  (1973)  have  demonstrated  that  also  FM— CW  radars  are 
capable  of  detecting  clear-air  convective  patterns. 


3.2.2  Radar  detection  of  stable  regimes 

In  stable  air  the  turbulence  and  also  the  refractive  index  fluctuations  are  suppressed  by 
negative  buoyancy  forces.  In  this  stable  regime  the  scatter  is  caused  by  refractive  index 
fluctuations  associated  with  turbulence  generated  by  mechanical  energy  extracted  from 
the  mean  flow.  Inside  regions  of  large  transfer  of  kinetic  energy  from  windshear  into 
turbulent  motion,  the  Cn2  can  according  to  Ottersten  (1969)  be  expressed  in  terms  of 
the  potential  refractive  index  gradient  d</>/dz,  the  gradient  of  horizontal  wind  du/dz  and 
the  flux  Richardson  number  Rj-  as 

Cn2  = a2e2/5(^)2(^-)  2|l-Rf]-‘  (3.1) 

from  Equation  (3.1)  we  sec  that  the  profiles  of  refractive  index  and  wind  may  create 
conditions  with  large  Cn2  and  enhanced  backscatter. 

3.2.3  Measurement  of  wind  field  using  radars 

The  wind  field  can  be  studied  using  Doppler  radars,  but  mainly  under  precipitation 
conditions  using  the  precipitation  particles  as  tracers. 

In  clear,  cloud-free  air,  a release  of  man-made  targets  is  usually  needed.  Such  targets 
might  be  thin  aluminium-coated  fiberglass  needles  cut  to  one-half  of  the  radar  wave- 
length employed. 

With  one  fixed  radar  only,  the  radial  wind  component  can  be  obtained.  From  a VAD 
scan  both  the  horizontal  wind  and  the  mean  profile  can  be  calculated. 

From  information  about  the  Doppler  spectral  variance,  information  on  the  turbulent 
kinetic  energy  may  be  derived  under  restricted  conditions.  The  variance  av2  of  the 
Doppler  spectrum  can  be  written  as 

°V2  = °S2+  aF2  + °B2  + <>T2 


where  the  individual  terms  are  the  variances  of  wind  shear,  os2;  the  fall  speed,  of2;  the 
finite  beamwidth,  oB2;  and  the  turbulent  kinetic  energy,  oT2.  To  separate  aT2  from 
this  equation  is  not  always  a simple  problem.  For  details  the  reader  is  referred  to 
Lhermitte  (1968)  and  Wilson  (1970). 


An  example  of  an  analysis  of  both  the  mean  wind  field  and  the  turbulence  field  is  shown 
in  Figure  3.4. 


WIND  SPEED,  ms  ' 


Figure  3.4  Mean  wind  field  and  turbulence  field  as  determined  from  * VI/)  Doppler 
information  in  snow  (from  Wilson.  1970) 

a)  temperature  sounding 

b)  mean  wind  profiles 

c)  turbulent  kinetic  energy  profiles 

d)  mechanical  production  rate 

The  velocity  resolution  is  better  than  0.25  ms  ' 


To  overcome  the  restrictive  time  stationarity  and  homogeneity  assumptions  which  have 
to  he  made  in  V AD,  two  or  three  radars  can  he  used  (COPLAN  scanning). 

An  example  of  the  potential  of  the  COPLAN  scanning  in  obtaining  two-dimensional 
wind  is  shown  in  Figure  3.5.  This  technique  is  powerful  in  examining  the  detailed 
structure  of  the  windficld  in  convective  storms. 

It  can  he  concluded  that  Doppler  radar  is  a powerful  tool  in  the  study  of  atmospheric 
motion. 
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hgure  J. 5 / wo-dimensional  vector  eddy  flow  field  obtained  from  the  radial  velocity 
fields  measured  by  two  Doppler  radars  operated  in  the  COPl.AN  scanning 
mode 

(Data  taken  by  the  radar  group  at  the  Wave  Propagation  Laboratory  in  Boulder.) 


3.3  Beyond  the  horizon  EM  scatter 

3.3.1  Determination  of  the  refractive  index  spectrum  >l'n(K.) 

Of  the  techniques  discussed  so  far,  forward  scatter  experiments  are  the  only  group 
which  has  been  successfully  used  to  deduce  the  three-dimensional  refractive  index 
spectrum. 

Ojessing  ( 1969)  has  discussed  the  sensitivity  of  different  forward  scatter  methods  to  the 
form  of  the  refractive  index  spectrum,  llis  results  are  shown  in  Figure  3.6,  where  he 
concludes  that  two-frequency  and  vertical  bcamswing  experiments  arc  the  most  sensitive 
to  the  form  of  the  refractive  index  spectrum. 

In  a later  paper  (Ojcssing,  1973)  hr  has  tried  to  interpret  the  variability  in  the  form  of 
the  measured  refractive  index  spectra  assuming  atmospheric  stratifications.  The  steeper 
spectra  may  be  explained  by  this  assumption,  llis  calculations  are  based  on  the 
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assumption  that  the  inertial  sub 
range  is  independent  of  the  layer 
thickness,  which  may  not  be  the 
rase  (sec  Metcalf  and  Atlas,  1973). 

Any  degree  of  anisotropy  and 
inhomogeneity  in  the  refractive 
index  field  will  complicate  the 
interpretation  of  the  data.  This 
can  to  a certain  extent  he  taken 
care  of  by  carefully  controlled 
synchronous  beam  swinging  in 
azimuth  and  in  elevation  of  both 
transmitter  and  receiver  beams 
(CJ jessing  and  Borresen.  1967). 


Hgure  3.6  I heoretical  relationship  between 
measured  quantity  ami  the  slope  n 
of  the  refractive- index  spectrum 
(^K)  ~ A,'  "j 

Length  of  the  arrows  shows  the  effect  of 
radiation  on  the  n deductions  (from  Cljess 
mg,  1969) 


3.3.2  Measurements  of  wind  by  forward  s»  alter  techniques 

When  discussing  wind  measurements  using  forward  scatter  Doppler,  one  must  distinguish 
between  *hc  post-detection  intensity  fluctuation  spectrum  S(F),  which  is  a measure  of 
the  Doppler  spread,  and  the  pre  detection  spectrum  S(f)  which  includes  the  Doppler 
shift  (Atlas.  1 969). 

I he  Doppler  spread  is  affected  by  a combination  of  the  cross-path  component  of  the 
wind  and  the  spatial  variation  of  the  wind  within  the  scattering  volume. 

II  we  have  access  to  the  absolute  Doppler  frequency  at  the  receiver,  we  can  obtain 
information  about  the  cross-path  wind.  Doppler  shift  (i  e the  mean  cross-path  wind)  can 
in  principle  be  measured  with  rake  delay  techniques  (Birkemcier.  1972). 

Recently  Olsen  and  I.ammcrs  (1973)  have  clearly  demonstrated  the  capability  of  a 
troposcatter  system  for  measuring  wind  velocity.  Some  of  their  results  from  a 500  km 
long  scatter  link  are  shown  in  Figure  3.7. 

For  further  details  regarding  the  use  of  F.M  forward  scatter  as  a remote  sensing 
technique  the  reader  is  referred  to  Cox  (1969). 


21 


TIME.  IN  HOURS  AFTER  0000  GMT.  10  JUNE  1970 


Figure  3.7  Comparison  of  scatter  and  radiosonde  measurements  of  cross-path  wind 
speed  at  6.9  and  12  km  using  rake  technique  (from  Olsen  and  l.ammers,  197.1) 
Scatter  measurements  at  10  km  off  path: 

o windward  side,  * leeward  side,  + radiosonde  measurements 


3.4  Probing  the  atmosphere  with  acoustic  sounders 

The  backscattered  acoustic  signal  is  determined  by  the  spectrum  of  the  acoustic 
refractive  index,  i e by  the  spectral  density  at  the  wavenumber  corresponding  to  one  half 
acoustic  wavelength.  This  is  directly  analogous  with  the  radar  case. 

For  pure  backscatter,  only  temperature  fluctuations  will  contribute,  while  in  the  case  of 
forward  scatter  (bistatic  mode)  both  the  wind  and  temperature  spectra  will  contribute 
to  the  scattering. 

Because  of  the  strong  influence  on  acoustic  waves  of  air  temperature,  acoustic  sounding 
technique  is  a powerful  tool  for  studying  atmospheric  structure  in  which  temperature 
fluctuations  can  be  used  as  tracers. 

As  in  the  case  of  radars,  two  different  regimes  of  atmospheric  turbulence  are  distin- 
guished; the  convection  regime  and  the  stable  regime. 

3.4.1  Detection  of  clear  air  convection 

Acoustic  echo-sounding  is  special- 
ly well  suited  for  studies  of  ther- 
mal activities  in  the  troposphere. 
As  an  example,  a typical  recording 
obtained  during  unstable  condi- 
tions is  shown  in  Figure  3.8.  The 
dark  portions  of  this  recording 
correspond  to  thermal  plumes. 
The  different  grey  tones  on  the 
facsimile  recordings  can  be  in- 
terpreted as  C-j-  values,  which 
again  are  a measurement  of  the 
temperature  fluctuations.  Rather 
good  agreement  (within  50%)  has 
been  found  between  sounder 
measurements  of  C.j  and  C.j.  cal- 
culated from  measurements  of 
temperature  differences  over  a dis- 
tance equal  to  one-half  acoustic 


IHCSI  ION 

Figure  3. S Facsimile  recording  of  backscatter 
acoustic  intensity  from  convective 
plumes  (from  Kjelaas  and  Ochs,  1974) 


wavelength  on  .1  nearby  meteorological  tower  (Hall,  1973).  Acoustic  sounding 
techniques  are  now  at  the  stage  were  can  In-  measured  as  a function  of  height  and 
comparison  can  be  made  with  different  existing  theories  for  convection. 


Detection  of  temperature  inversions 

l'wo  types  of  temperature  inversion  or  stable  layers  may  be  detected  by  acoustic- 
sounding  technique.  On  clear  nights,  thermal  radiation  from  the  earth’s  surface  results  in 
a lowering  of  the  ground  temperature.  Turbulence  or  molecular  conduction  causes  a 
heat  loss  to  the  earth  and  a layer  of  cold  air  forms  just  above  the  earth’s  surface.  As  the 
heat  conduction  losses  increase  during  the  night,  the  depth  of  this  cold  layer  increases. 
In  a case  like  this  we  would  not  expect  a random  thermal  structure  and  the  acoustic- 
sounder  would  only  detect  the  temperature  gradient  separating  the  cold  air  from  the 
warmer  air  above  (Hall,  1973).  A second  type  of  temperature  inversion  is  caused  by 
subsidence. 

A typical  facsimile  display  illustrating  the  potential  of  the  technique  under  stable 
conditions  is  shown  in  Figure  3.9.  The  dark  portions  of  the  picture  are  regions  from 
which  echoes  were  received.  In  this  figure  different  types  of  atmospheric  structures  are 
illustrated,  lhc  lowest  dark  portion  of  the  records  is  the  nocturnal  inversion.  At  about 
0600  hours  this  radiation  inversion  starts  lifting  and  finally  at  about  1000  hours  free 
convection  is  more  or  less  developed.  Above  this  inversion  we  observe  a wavy  structure 
which  is  probably  caused  by  wind  shear  induced  instabilities  or  Kelvin  Helinholtz 
billows.  Hook  et  al  (1972)  have  reported  .ui  excellent  agreement  between  the  observed 
K-  H waves  or  gravity  waves,  and  the  corresponding  pressure  fluctuation  observed  on  a 
microbarograph.  Kjelaas  et  al  (197-la)  have  demonstrated  the  capability  of  an  array  of 
acoustic  sounders  in  determining  phase  velocities  of  gravity  waves. 
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Figure  3.9  Facsimile  recording  o)  backseat  ter  acoustic  intensity  during  conditions  of 
nocturnal  inversion  and  breaking  up  of  the  same 

(Data  taken  by  the  Acoustic  (iroup  of  the  Wave  Propagation  Laboratories,  Boulder.) 


Under  stable  conditions  the  agreement  between  the  Cr  derived  from  acoustic  sounding 
measurements  and  C(.  calculated  from  temperature  measurements  is  not  as  good  as 
when  the  atmosphere  is  convectively  unstable.  Ibis  may  be  due  to  reflection  or 
enhancement  of  the  returned  signal  because  of  focusing  effects. 


Measurements  of  wind  velocity  using  .uoiistic  sounding  technique 

As  in  the  case  « >t  KM  propagation,  the  Doppler  frequency  shift  can  he  used  to  determine 
the  wind  component  along  the  dn.(ks  k0),  whii  It  is  the  bisector  of  the  angle  formed  l>\ 
the  itu  ident  heam  and  the  scattered  (ream,  for  acoustic  waxes,  refraction  of  the  beam 
may  be  a problem  and  can  lead  to  misinterpretations.  An  analysis  of  the  influence  of 
iH'am  bending  is  therefore  important  for  every  experimental  configuration  of  acoustic 
soundings  (Derail  and  Clifford,  I‘*7L’).  A typical  Doppler-derived  vertical  velocity  pat 
tem  for  a convective  plume  is  shown  in  f igure  A 10.  On  the  left  hand  side  of  t fie  figure 
is  the  corresponding  b.u  kscatter  intensity. 
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hgurej.K)  Acoustic  echoes  (left)  and  Doppler  derived  vertical  velocity  field  during 
thermal  plume  at.  ' -tVy  (from  Horan  et  al,  1971) 

C.ontour  interval  is  0.4  ins  '.  regions  of  vertical  velocity  >0.8  ms  1 are  shaded;  arrows 
indicate  direction  ot  vcrtu.il  i iMiiponcnt. 


figure  .1.11  shows  an  example  of  the  total  wind  vector  determined  b\  employing  three 
inclined  acoustn  echo  sounders.  1 fie  agreement  xxith  the  radiosonde  derived  wind 
velocity  is  convitu  mg;  the  teclmi<|ue  has  great  possibilities  as  a wind  sensor  in  the  lower 
troposphere. 
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Figure  4.11  Comparison  between  wind  profile  measured  with  radiosonde  and  the  one 
derived  from  acoustic  Doppler  measurement  ( from  Beran  et  al,  1974) 

3.4.4  Flux  measurements 

Using  similarity  theory,  which  seems  to  model  free  convection  well,  heat  flux  can  be 
related  to  CrJ;  accordingly  heat  flux  can  be  gained  from  the  backseat tered  acoustic 
energy. 

Using  intersecting  beams  from  three  sounders,  the  wind  vector  can  be  measured.  When 
velocity  changes  from  pulse  to  pulse  are  computed,  thus  yielding  u'w'  and  v'w',  a 
measure  of  the  vertical  momentum  flux  or  stress  is  derived  (Hall,  private  communica- 
tion, 1973). 

Further,  the  dissipation  may  be  measured  by  analyzing  pulse  to  pulse  velocity  dif- 
ferences, the  second  moment  of  the  Doppler  spectra,  or  alternatively  from  bistatic 
scattering  measurements,  where  Cv  3 contributes  to  the  acoustic  scatter. 

3.5  Summary  of  experimental  techniques 

Table  3.1  summarizes  the  relative  merits  of  the  techniques  discussed  in  the  previous 
sections.  Again  we  want  to  emphasize  that  only  the  techniques  where  the  received  signal 
is  determined  by  refractive  index  inhomogenities  are  discussed.  Thus  factors  such  as 
absorption  and  molecular  scattering  where  the  signal  reacts  with  individual  molecules  are 
excluded. 

From  Table  3.1  we  sec  that  no  single  technique  provides  information  about  all  the  listed 
parameters  for  describing  the  atmospheric  structure.  If  more  than  one  technique  is  used 
simultaneously,  we  shall  have  to  remember  that  the  various  techniques  may  be  sensitive 


to  iliffcmit  portions  of  the  refractive  index  wavenumber  spe.  (rum.  ic  different  s.  ale 
si/e  ot  the  refractive  index  inhomogcnilics,  viz 

radar  and  acoustic  ~ I 10  cm  (X/2) 

optical  line-of-sight  ~ cm  (first  Fresnel  zone) 

microwave  line-of-sight  ~ a few  in  (first  Fresnel  zone) 

forward  scatter  ~l  20  tn  (X/2  sin  0/2) 

In  addition  to  being  sensitive  to  different  s».de  sizes,  the  variables  measured  max  he 
averaged  over  different  distances  and  volumes. 

For  line-of-sight  the  average  volume  is  determined  hy  the  length  of  the  path  and  some 
weighting  along  the  path,  and  In  the  I team  width. 

In  the  case  ol  radar  and  acoustic  backs,  altering  devices,  the  volume  over  which  the 
variable  is  averaged  is  determined  hv  the  heamwidth  and  bv  the  length  of  the  trans 
milted  pulse. 

For  forward  scatter  the  volume  is  given  In  the  beam  width,  scattering  angle  and  distance 
between  transmitter  and  receiver. 

We  shall  now  briefly  discuss  Fable  3.  I. 

3.5.1  Wind 

From  Table  3.1  we  see  that  mean  wind  ;uid  wind  fluctuation  can  be  measured  in  several 
different  ways,  although  in  the  absence  of  refractive  index  fluctuations  of  sufficient 
intensity  a tracer  such  as  snow,  rain  or  chaff  must  be  employed  in  the  radar  hackscatter 
case.  1 his  has  been  the  case  for  Doppler  radar  studies,  where  a two-dimensional  coverage 
ot  the  wind  field  has  been  obtained.  This  is  not  the  case  with  forward  scatter  when  the 
sensitivity  is  (sin  0/2)'"  ' better  than  for  backs,  alter. 

A line  average  of  the  transverse  wind  can  be  obtained  from  both  microwave  and  optical 
line-of-sight  propagation.  The  two- dimensional  horizontal  wind  has  also  been  success- 
fully measured  with  optical  line  of  sight  propagation  methods.  The  vertical  profile  of  the 
three-dimensional  wind  can  be  obtained  using  aconsti.  sounding  Doppler  techniques. 
From  acoustic  sounding  some  knowledge  about  1 , the  turbulence  structure  constant, 
can  be  gained.  As  acoustic  technique  also  gives  the  three  dimensional  wind,  momentum 
flux  can  also  lie  derived.  Profiles  of  transverse  wind  can  also  be  measured  from  F.M 
forward  scatter  techniques,  providing  absolute  Doppler  can  be  measured. 

from  optical  line-of-sight  propagation  techniques,  horizontal  divergence  has  been 
successfully  derived. 

3.5.2  Temperature 

Because  of  the  strong  dependence  of  humiditv  rather  than  of  temperature  on  the  radio 
refractive  index,  available  techniques  for  temperature  probing  are  limited  to  optical  and 
acoustic  waves.  A line  average  of  the  mean  temperature  as  well  as  of  the  intensitv  of  the 
temperature  fluctuation,  (l..3,  can  be  derived  from  optical  line  of  sight  propagation 
measurements.  In  some  cases  also  information  about  the  shape  of  the  temperature 
spectrum  can  be  obtained. 

Acoustic  sounding  techniques  can  also  give  a line  profile  of  the  C'. , 3 averaged  over  some 
volume  (pulse  volume).  During  adiabatic  conditions  the  vert  it  al  heat  flux  can  be  derived 
from  the  C.j 3 and  the  vertical  wind  velocitv. 
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3.5.3  Humidity 

Measurements  of  humidity  are  difficult  to  perform  because  the  optical  and  acoustic- 
refractive  indexes  are  insensitive  to  humidity  and  the  microwaves  are  influenced  both  by 
temperature  and  humidity.  By  assuming  the  temperature  more  or  less  constant  over  a 
path,  microwave  line-of-sight  technique  has  been  used  to  measure  a line  average  of  the 
water  vapour  pressure. 


3.5.4  Refractive  index 

As  discussed  in  section  2.1,  the  refractive  indexes  for  optical,  micro,  and  acoustic  waves 
are  mostly  dependent  on  temperature,  humidity  and  temperature,  and  temperature  and 
wind  velocity,  respectively.  Microwave  line-of-sight  will  give  a line  average  of  both  the 
mean  and  the  fluctuations  Cn2  of  the  radio  refractive  index,  and  in  some  cases 
information  about  the  form  of  the  refractive  index  spectrum. 

For  backscatter,  both  microwave  and  acoustic,  the  intensity  of  the  refractive  index 
fluctuations  Cn2  can  be  measured  as  a function  of  height.  Microwave  forward  scatter  has 
been  used  to  derive  the  three-dimensional  refractive  index  spectrum. 

As  noted  earlier,  different  techniques  are  sensitive  to  different  scale  sizes,  averaging 
distances  and  volumes  and  in  most  cases  arc  suitable  for  only  limited  regions  of  the 
atmosphere.  In  Table  3.2  we  have  divided  the  atmosphere  into  three  regions  and  listed 
the  techniques  which  may  be  used  within  each  region: 

Line-of-sight  — mainly  suitable  for  measurements  within  a few  meters  above  the  ground, 
although  paths  from  sea-level  to  mountain-top  have  been  used,  and  earth— satellite 
measurements  are  becoming  available. 

Radars  — different  radars  have  been  used  from  a few  meters  above  ground  to  many 
kilometers,  both  vertical  and  slant. 

Acoustic  sounding  — greatest  potentials  from  a few  meters  up  to  a few  kilometers. 

Forward  scatter  — since  generally  beyond  horizon,  usable  region  from  a few  kilometers 
above  ground  to  many  kilometer*. 
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Table  3.2  Height  intervals  over  which  different  remote  sensing  technioues  max  he  used 
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CONCLUSION 
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It  can  be  concluded  that  t art1  fully  designed  wave  propagation  and  m alter  experiments 
can  provide  information  about  the  atmospheric  strut  turr  otherwise  not  readily  obtain 
able.  Some  of  the  techniques  discussed  in  this  paper  are  in  a formative  rather  than  a 
productive  stage.  This  is  particularly  the  case  with  forward  scatter  and  line  of  sight 
experiments,  although  optical  line -of  sight  seems  to  be  a very  useful  remote  wind  sensor. 
Acoustic  sounding  techniques  are  in  a stage  of  rapid  development  and  it  has  been 
demonstrated  that  this  technique  has  very  great  potentials  in  sensing  both  wind  and  the 
thermal  structure  of  the  atmosphere.  High  resolution  radars,  such  as  the  KM  CW  radar, 
have  revealed  the  existence  of  thin, undulating  layers  and  waves  in  the  atmosphere  which 
have  not  been  known  before.  These  waves  (K  II  waves)  may  be  of  importance  for  a 
further  understanding  of  the  transport  of  momentum  in  the  houndarv  layer. 

The  high  power  pulse  radars  contribute  to  the  understanding  of  t leal  ait  ei  ho  (('A  T)  and 
the  study  of  gravity  waves.  Multistation  pulse  Doppler  radars  scanning  in  a plane 
(COPLAN)  are  used  to  study  the  horizontal  wind  field  within  thunderstorms,  snow 
storms  and  also  in  the  clear  .ur  when  employing  chaff  (artificial  tracer).  These  expert 
ments  add  new  and  useful  information  to  the  know  ledge  about  the  cell  strut  lure  ol  sin  h 
storms. 
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ABSTRACT 

A wind-measuring  system  using  three  separate  optical  systems,  each  consisting  of  a 3 m\V  He-Nc  laser 
and  two  photodiode  receivers,  forming  an  equilateral  triangle  300  m on  a side,  has  successfully  measured 
the  divergence  over  the  area  of  the  triangle,  and  the  space-averaged  horizontal  wind  vei  tor.  Good  correla- 
tion was  found  between  the  How  into  the  triangle  and  occurrences  of  thermal  plumes  seen  by  an  adjacent 
acoustic  sounder.  The  How  into  the  triangle  was  proportional  to  the  vertical  velocity.  I hiring  large  con- 
vective activity,  there  was  a certain  periodicity  in  the  occurrence  of  plumes. 


1.  Introduction 

Reliable  measurements  of  the  divergence  of  the 
horizontal  wind  are  difficult  to  achieve.  Many  attempts 
have  been  made  to  develop  methods  based  on  point 
measurements  of  wind.  On  the  synoptic  scales, 
Bellamy  (1949)  developed  a technique  that  used  a 
triangle  of  wind  reports  to  calculate  the  divergence  for 
the  enclosed  area;  Kndlich  and  Clark.  (1963)  calculated 
divergence  by  linear  interpolation  between  wind  com- 
ponents at  the  vertices  of  the  triangle  of  observation. 
Eddy  (1964)  improved  Bellamy’s  triangle  technique 
bv  using  a computer  model.  He  also  applied  a correla- 
tion technique  to  establish  the  space  correlation  and 
time  persistence. 

On  smaller  scales  (i.e.,  those  scales  comparable  with 
plumes  and  on  up  to  a few  hundred  meters),  no  reliable 
techniques  for  measuring  the  horizontal  divergence 
exist.  Many  laboratory  experiments  have  formed  the 
basis  for  different  theories  for  entrainment  into  plumes. 
Morton  ft  al.  (1956)  have  described  a model  assuming 
that  the  rate  of  entrainment  into  the  plumes  is  propor- 
tional to  the  vertical  velocity  inside  the  plume.  Telford 
(1966)  has  developed  a model  assuming  that  the  entrain- 
ment should  not  be  related  to  the  mean  velocity  but  to 
the  local  level  of  turbulence. 

Due  to  lack  of  suitable  wind -measuring  techniques 
these  theories  have  been  difficult  to  check  in  the  atmo- 
sphere. Recently,  remote  sensing  of  space -averaged 
wind  has  progressed,  making  it  possible  to  calculate  the 

1 The  paper  was  prepared  while  the  senior  author  was  a visiting 
scientist  at  the  National  Oceanic  and  Atmospheric  Administra- 
tion, Environmental  Research  Laboratories,  Boulder,  under 
sponsorship  of  the  Royal  Norwegian  Council  for  Scientific  and 
Industrial  Research. 


divergence  for  a small  enclosed  area.  Lawrence  el  al. 
(1972)  have  described  an  optical  wind  sensor  that  uses 
the  naturally  occurring  density  fluctuations  to  obtain 
a spatially  averaged  wind.  By  using  a correlation 
technique  at  the  receiver,  they  are  able  to  extract  a 
path  averaged  transverse  wind  speed. 

This  paper  discusses  field  measurements  of  divergence 
;ind  wind  velocity  based  on  an  extension  of  the  optical 
technique.  A He  N'e  laser  triangle  path  (300  m on  a 
side)  was  set  up  near  a 150  m tower.  The  results  are 
compared  with  wind  velocity  measured  at  five  levels 
on  the  tower  and  with  the  occurrence  of  plumes  as  seen 
by  an  acoustic  sounder  located  inside  the  triangle. 

2.  The  optical  measuring  equipment 

The  optical  technique  used  for  wind  measurement  is 
described  fully  in  Lawrence  el  al.  (19/2);  however,  a 
brief  discussion  is  included  here  for  completeness,  since 
it  is  a recent  development.  The  technique  may  be  used 
over  ranges  from  3lX)  m to  more  than  15  km,  so  that  it 
is  particularly  useful  for  obtaining  spatially  averaged 
wind  measurements.  In  addition,  it  is  effective  for 
measuring  low  mean  wind  flow  in  the  presence  of  large 
fluctuating  components. 

The  technique  uses  the  drifting  scintillation  pattern 
arising  from  wind-transported  refractive-index  irregu- 
larities in  the  atmosphere.  This  drifting  scintillation 
pattern  is  the  cause  of  the  twinkling,  or  scintillation,  of 
starlight.  With  a very  small,  bright  light  source,  such 
as  a laser  shining  through  the  atmosphere,  the  moving 
scintillation  pattern  can  be  seen  clearly  on  a screen  at 
night.  The  sizes  of  the  irregular  patches  of  light  in  the 
pattern  are  determined  by  the  diffraction  of  the  incident 
beam;  they  result  from  refractive-index  irregularities 
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(almost  entirely  the  result  of  temperature  differences) 
at  various  |H>sitiotis  along  the  propagation  path.  At 
each  position,  the  diffraction  process  selects  a most 
effective  stale  of  irregularity,  or  eddy,  whose  size  is 
sufficient  to  cause  a half  wavelength  difference  in  path 
length  to  the  receiver  for  two  extreme  rays  passing 
through  its  edges.  The  diffraction  effects  of  this  par- 
ticular eddy  dominate  over  those  of  other  eddies  at  the 
same  location. 

It’  a diverging  laser  beam  is  effectively  a point  source, 
geometric  magnification  increases  the  si/e  of  the  received 
eddy  image  by  the  factor  /.  s,  where  z is  the  [Position 
along  a beam  of  length  The  image  drift  rate  at  the 
receiver  caused  by  the  local  eddy  drift  velocity,  de- 
pending on  the  wind  at  position  2,  is  also  increased  by 
the  same  factor.  Therefore,  the  receiver!  optical  distur 
bance  is  larger  and  moves  faster  than  the  eddy  that 
produced  it.  Since  different  portions  of  the  path 
contribute  to  the  total  scintillation  pattern,  that  pattern 
will  be  continuously  evolving  because  of  the  relative 
speeds  of  the  different  sized  features  it  contains. 

To  measure  the  average  component  of  the  wind 
across  an  optical  path,  we  look  at  fluctuations  in  light 
intensity  in  two  small  apertures  spaced  along  the  direc- 
tion of  movement  of  the  pattern,  as  shown  in  Fig.  1. 
A function  like  that  shown  in  Fig.  2 results  if  the  nor- 
malized covariance  of  the  logarithms  of  these  signals  is 
computed.  For  the  conditions  of  a uniform  wind  at 
right  angles  to  the  optical  path  and  parallel  to  the 
separation  of  the  detectors,  and  for  uniform  refractive- 
index  structure  constant  (CV)  along  the  path,  some 
measure  of  the  time  scale  of  the  function,  e.g.,  the  time 
delay  to  the  peak  or  the  slope  of  the  function  at  zero 
delay,  could  be  used  to  derive  wind  speed.  For  the  real 
case,  with  winds  of  arbitrary  direction  and  speed  along 
the  path,  and  with  a non  uniform  distribution  of 
along  the  path,  this  function  changes  not  only  in  time 
scale,  but  also  in  shape.  It  can  be  shown  (Lawrence  el  ol., 
1972),  for  uniform  CV  along  the  path,  that  the  slope  of 
the  function  at  zero  delay  has  a known  relationship  to 
the  average  horizontal  wind  across  the  path,  even 
though  the  wind  along  the  path  is  non  uniform  in 
direction  and  speed.  Thus,  we  measure  the  slope  at  zero 
delay  to  determine  the  average  wind.  The  weighting 
function  of  the  wind  measurement  for  various  receiver 
separations  is  shown  in  Fig.  ,1  as  a function  of 
d“g(\/.)  *,  where  p is  the  detector  spacing,  X the  laser 
wavelength,  and  l.  the  path  length.  For  the  divergence 
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Kio.  I.  The  general  arrangement  used  to  measure  the  average 
component  of  the  wind,  across  a laser  beam,  parallel  to  the 
(usually  horiiontal)  spacing  of  two  detectors. 
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Kio  2 A schematic  diagram  of  the  normalized  covariance 
function,  showing  the  delay  to  the  |ieak  and  the  slo|>e  at  zero 
time  lag  While  the  former  is  frequently  used  to  measure  pattern 
drift  velocities,  the  latter  is  more  suitable  in  the  optical  case 
where  the  pattern  decays  rapidly. 

measurement  reported  here,  /S— 0.46.  This  separation, 
and  the  use  of  apertures  t>  mm  in  diameter,  gave  a 
nearly  symmetrical  weighting  function  around  the 
center  of  the  path. 

There  are  several  limitations  to  the  present  optical 
system,  but  we  no  not  believe  that  they  affected  the 
wind  measurements  reported  here.  The  wind  weighting 
function  is  strictly  correct  only  when  the  wind  measure 
nient  is  averaged  over  enough  time  so  that  CV  is 
statistically  uniform  over  the  path.  ’This  can  be  satisfied 
for  a horizontal  path,  but  if,  for  example,  a slant  path 
is  used  (V  will  tend  to  be  higher  at  the  lower  end  of  the 
luth,  and  the  wind  in  this  portion  will  be  weighted 
more  heavily.  Random  variations  of  (. V will  distort  the 
measured  average  wind  sjieed  only  if  they  are  correlated 
with  the  variations  in  wind  speed.  A reasonable  sup- 
(Hisition  is  that  turbulence  is  strongest  at  locations  of 
maximum  wind  shear,  not  at  locations  of  maximum 
wind  speed.  If  so,  the  turbulence  tends  to  lie  uncorre- 
lated with  the  wind,  and  the  average  obtained  by  the 
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Kio.  ,t.  The  relative  weights  of  the  different  portions  of  the 
path  in  determining  the  optically  measured  wind.  The  parameter 
jS-atXZ,)'*  is  the  separation  of  the  sensors.  These  curves  are 
calculated  for  point  sensors;  the  Unite  area  of  a real  sensor  will 
tend  to  remove  the  negative  weights  near  the  receiver. 
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Flo.  4.  A comparison  of  the  average  wind  speed  (I’d)  measured  optically,  with  the  average 
of  the  readings  of  six  Gill  propellor  anemometers  (!'.)  uniformly  spaced  along  the  optical 
path.  Both  the  anemometers  and  the  optical  system  measure  the  component  of  wind  that  is 
horizontal  and  at  right  angles  to  the  path.  The  path  is  300  in  long,  oriented  north  south,  and 
3 m above  the  ground. 


optical  system  is  unbiased.  A second  source  of  error 
may  occur  when  there  is  a combination  of  high 
refractive-index  fluctuations  and  long  path.  Under  these 
conditions,  a so-called  saturation  effect  occurs,  in  which 
the  optical  scintillation  no  longer  increases  with  C.1. 
When  this  condition  occurs,  the  spherical  propagation 
equations  upon  which  the  optical  wind  measurement  is 
based  are  no  longer  correct.  The  result  is  that  the  optical 
measurement  will  underestimate  the  true  wind,  and  in 
severe  saturation  this  may  be  by  a factor  of  2 or  more. 
Saturation  may  be  easily  detected  by  the  optical  system 
and  did  not  occur  on  the  300  m paths  used  in  our 
divergence  experiment. 

Before  undertaking  the  divergence  experiment,  we 
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Fio.  5.  The  layout  of  the  laser  triangle,  130  in  meteorological 
tower,  and  the  acoustic  sounder. 


compared  the  average  wind  speed  reading  of  one  of  the 
three  optical  systems,  over  the  same  height  and  path 
length  that  was  later  used,  with  the  average  of  the 
readings  of  six  Gill  propellor  anemometers  uniformly 
spaced  along  the  optical  path.  The  anemometers  were 
oriented  to  read  the  component  of  wind  that  was 
horizontal  and  at  right  angles  to  the  optical  path.  The 
result  for  a 24  hr  period  is  shown  in  Fig.  4 and  appears 
to  otter  full  proof  of  satisfactory  performance. 

3.  Measurement  and  data  reduction 

The  measurements  were  made  at  Haswell,  Colo., 
during  a two-week  period  in  the  beginning  of  August 
1972.  The  Haswell  area  has  a typical  high-plain 
meteorological  regime  with  strong  radiative  cooling 
during  the  night  and  strong  solar  heating  and  convective 
activity  during  the  day. 

Three  separate  optical  systems,  each  consisting  of  a 
3-mW  He-N'e  laser,  two  photodiode  receivers,  and  a 
small  specialized  computer  providing  a real-time  read- 
out of  average  wind  speed,  were  used  in  the  experiment. 
The  optical  paths  formed  an  equilateral  triangle,  300  m 
on  a side  and  2.3  m above  the  ground.  Fig.  5 shows  the 
location  of  the  laser  triangle  relative  to  the  meteoro- 
logical tower  and  the  acoustic  sounder.  The  terrain  near 
the  triangle  is  flat  and  free  from  any  obstacles  for  a 
minimum  of  3 km  in  any  direction.  The  area  is  covered 
with  clumps  of  buffalo  grass  (^15  cm  in  height).  Since 
each  optical  path  was  arranged  to  read  an  outflowing 
wind  as  positive,  the  sum  of  the  readings  of  all  three 
paths  represents  the  average  horizontal  wind  sjwed  out 
of  the  triangle.  The  individual  wind  tradings  for  each 
optical  path,  smoothed  with  a Id  sec  running  average, 
as  well  as  the  scalar  sum  of  all  three,  were  recorded 
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Locot  Tim« 

Kio  6.  Ten-second  Average  divergence  recordings  for  typical  afternoon  and  night  periods. 


simultaneously  on  strip  chart  recorders  and  on  a digital 
tape  recorder. 

The  horizontal  wind  divergence  is  given  bv 

/.(t'o+t'uo+t’sso) 

V»V- , 

.4 

where  ro,  riW,  vf<1)  are  the  wind  components  transverse 
to  the  east-west,  30°  210°  and  150°  -330°  triangle  sides, 
respectively ; l.  is  the  side  of  the  triangle;  and  .4  the 
area  of  the  triangle. 

Typical  recordings  of  divergence  during  periods  of 
large  convective  activity  and  stable  conditions  are 
shown  in  Fig.  6.  Convergence  of  10  1 sec  1 over  periods 
of  minutes  is  common  during  daytime.  During  the 
night,  values  of  5X10-*  sec  1 are  normal.  This  value 
may  be  considered  the  upper  limit  of  the  noise  level  of 
the  divergence  measurement.  In  the  comparison  be- 
tween divergence  and  vertical  velocity,  a 100-sec  run- 
ning average  was  used.  The  space-time  average  wind 
velocity  was  computed  by  using  the  wind  components 
measured  over  two  of  the  triangle  sides.  This,  of  course, 
gives  a redundant  estimate  of  the  velocity  as  we  have 
three  pairs  of  legs  on  a triangle. 

Using,  for  example,  the  east  west  and  30°  210°  pair, 
the  horizontal  space  average  wind  speed  T,  can  be 
expressed  in  terms  of  the  two  components  r0  and  t>1M  as 


and  the  direction  a from  which  the  wind  is  blowing  as 

/ t>0v3  \ 

a-tan-'f )+180°. 

'2  V uoftv 


Using  the  other  pairs  of  the  triangle,  similar  expressions 
are  obtained. 

The  meteorological  tower  was  instrumented  to  mea- 
sure vector  wind  by  the  use  of  bivanes  at  five  fixed 
levels:  31,  62,  93,  124  and  140  m.  The  acoustic  sounder 
inside  the  triangle  was  pointing  vertically,  and  its 
records  were  compared  with  the  measured  inflow  into 
the  triangle.  All  data,  except  the  acoustic  sounder  data, 
were  sampled  once  per  second  and  recorded  on  digital 
ta[>es.  The  acoustic  sounder  data  were  recorded  (both) 
on  facsimile  and  analog  tapes. 

4.  Results 

tj.  Comparison  of  divergence  and  the  acoustic  sounder 

Considerable  evidence  now  exists  that  plume  like 
structure  can  be  seen  by  an  acoustic  sounder  (McAllister 
rt  al.,  196*7;  Heran,  1971).  Fig.  7 shows  an  acoustic 
sounder  record  and  a divergence  for  predominately 
large  convective  activity.  The  dark  portions  on  the 
sounder  record  correspond  to  thermal  plumes  and  their 
height  indicates  the  height  of  the  plumes,  while  the 
white  “gaps"  are  interpreted  as  area  between  plumes. 
The  solid  line  is  the  divergence  measured  by  the  optical 
system.  Convergence  is  upward  on  the  figure. 

The  correspondence  between  the  measured  conver- 
gence and  the  plumes  as  seen  by  the  acoustic  sounder  is 
mostly  good,  although  some  discrepancies  can  be  seen 
especially  after  1440  local  time.  Note  that  while  the 
divergence  is  measured  at  a height  of  about  2 m above 
ground,  the  lowest  height  seen  by  the  acoustic  sounder 
is  about  50  m,  making  an  absolute  comparison  between 
the  two  systems  difficult. 
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Fig.  7.  Comparison  1*1  ween  the  plumes  seen  by  the  acoustic  sounder  and  the  flow  into 
the  triangle.  Large  plumes  are  associated  with  large  convergence.  The  solid  line  is  the  diver- 
gence measured  by  the  optical  system. 


b.  Correlation  of  divergence  anti  vertical : find 

The  vertical  velocity  is  related  to  the  divergence  field 
through  the  equation  of  continuity  as 

chr 

V»-V- . 

dt 

Fig.  8 shows  a typical  set  of  cross-correlation  curves 
between  the  divergence  and  the  vertical  wind  and 
gradient  on  the  tower  during  strong  convective  activity 
Because  the  tower  is  outside  the  triangle,  the  peak  in 
the  cross  correlations  does  not  occur  at  zero  lag  but  is 


Log  in  Seconds 


Fio.  8.  Cross-correlation  curves  between  divergence  and  vertical 
wind,  and  divergence  and  vertical  wind  gradient  for  the  different 
levels  on  the  tower.  Solid  lines  are  divergence  vs  vertical  velocity 
and  dotted  lines  divergence  vs  vertical  wind  gradient. 


shifted  in  time  corresponding  to  the  travel  time  between 
the  tower  and  the  center  of  the  triangle.  For  the  case 
chosen  the  wind  direction  was  along  a line  between  the 
center  of  the  triangle  and  the  tower.  The  negative  lag 
means  that  the  divergence  is  lagged  behind  the  vertical 
wind  on  the  tower. 

One  of  the  interesting  features  in  Fig.  8 is  the  differ- 
ence in  correlation  between  T»  V and  —dw/dt,  and 
between  Y»  ■ V and  w.  From  these  curves,  we  see  that  in 
the  case  of  plumes,  the  convergence  into  the  plumes  is 
more  determined  by  the  vertical  wind  than  by  the 
gradient  of  the  same.  These  results  seem  to  agree  with 
Morton  et  al.  (1956),  who  assumed  that  the  entrainment 
into  plumes  was  proportional  to  the  upward  velocity. 
It  might  be  argued  that  this  picture  might  have  been 
different  if  we  had  had  a measurement  of  the  vertical 
wind  at  the  same  height  as  the  divergence  measure- 
ment. Pue  to  the  100-sec  running  average,  however, 
we  are  looking  at  a large  population  of  plumes  which 
should  not  differ  much  from  5 to  30  m,  so  we  believe 
that  our  results  give  a reasonably  correct  picture. 

Another  interesting  feature  of  these  curves  is  the 
relative  time  shift  between  the  peaks  in  the  cross 
correlations  for  the  vertical  wind  at  the  different  levels 
indicating  that  the  plumes  or  convective  currents  are 
tilted.  The  tilt  is  in  the  downwind  direction.  This  corre- 
sponds to  what  Kaimal  and  Businger  (1970)  have 
reported  for  a single  plume.  The  plumes  are  passing  the 
four  upper  levels  on  the  tower  at  about  the  same  time 
but  some  20-30  sec  earlier  than  the  lowest  level  (30  m). 

In  all  the  cases  we  have  analyzed,  the  plumes  were 
tilted  in  the  lowest  60  m,  whereas  from  60  to  150  m 
almost  no  tilting  was  observed. 

The  third  feature  that  should  be  observed  in  Fig.  8 is 
the  periodicity  in  the  cross-correlation  curves  indicating 
that  there  is  a periodicity  in  the  development  of  con- 
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vective  currents  in  the  boundary  layer.  The  time 
between  the  two  cross-correlation  maxima,  which  might 
be  interpreted  as  the  period  between  a set  of  plumes, 
is  360  sec. 

Fanofsky  (1973)  has  recently  shown  that  the  plumes 
seem  to  move  with  a speed  close  to  the  local  wind.  In 
this  case  the  wind  was  about  1.5  n-  sec-1  at  the  height 
of  the  triangle,  giving  us  a separatvin  in  space  of  about 
540  m between  every  major  set  o;  plumes.  This  scale 
size  seems  to  agree  with  what  Aoierman  (1967)  has 
reported  for  convective  clouds.  She  found  scale  sizes 
of  400-900  m.  This  periodicity  was  only  found  during 
large  convective  activity.  During  the  morning  and  the 
evening  hours  correlation  curves  like  curve  B shown  in 
Fig.  9 are  typical.  The  magnitude  of  the  correlation 
between  TV  V and  w seems  to  change  very  little  during 
the  day,  but  the  plumes  are  released  more  randomly 
with  no  outstanding  size  contributing  to  the  cross 
correlation  during  the  morning  and  evening  hours. 

c.  Daily  variation  of  the  divergence 

As  shown  in  Fig.  6,  a great  difference  exists  between 
the  divergence  during  day  and  night.  In  the  morning 
and  evening,  we  found  that  the  flow  was  out  of  the 
triangle  about  60%  of  the  time.  During  heavier  con- 
vective activity  in  the  middle  of  the  day,  the  flow  was 
into  the  triangle  60%  of  the  time.  In  the  middle  of  the 
day  peak  values  in  convergence  (10-*  sec-1)  were  found 
over  periods  of  1-2  min.  The  standard  deviation  was 
normally  in  the  order  of  4X10™*  sec-1. 

d.  Space  average  of  icin,l  velocity 

As  shown  in  Fig.  4,  this  technique  measures  the 
transverse  space-averaged  wind  component.  By  com- 
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Fio.  9.  Cross-correlation  curves  between  divergence  and  vertical 
wind.  The  solid  line  is  for  high  convective  activity,  and  the 
dotted  line  for  low  convective  activity. 


Local  time 

Kio.  10.  A comparison  of  the  three  values  of  the  wind  vector 
determined  from  the  laser  triangle.  Wind  direction  and  speed  are 
shown  in  the  upper  and  lower  parts  of  the  figure,  respectively. 
A 20-min  running  average  has  been  applied  to  the  data  to  smooth 
out  the  effect  of  convective  activity. 

bining  two  components  as  shown  in  Eqs.  (2)  and  (3), 
we  get  the  total  horizontal  wind  vector.  Since  three 
components  are  available,  these  can  be  combined  to 
give  three  values  of  the  wind  vector.  By  applying  a 
20-min  running  average  to  the  data  to  smooth  out 
convective  activity  and  comparing  these  three  values, 
we  have  a measure  of  the  ability  of  the  system  to  obtain 
the  space-time  averaged  horizontal  wind  vector. 

Fig.  10  shows  plots  of  the  three  values  of  the  wind 
vector  for  nearly  6 hr,  with  direction  and  speed  shown 
in  the  upper  and  lower  parts  of  the  figure,  respectively. 
For  a wind  speed  higher  than  about  2 m sec-1,  we  see 
that  it  makes  little  difference  which  two  of  the  three 
components  are  used  to  calculate  the  horizontal  wind 
vector.  When  the  wind  speed  drops  below  about 
2 m sec-1  the  differences  increase,  very  likely  from 
several  causes.  Due  to  slight  differences  in  the  terrain 
along  the  triangle  legs,  irregular  wind  patterns  can 
persist  over  the  triangle  and  be  more  pronounced  during 
light  wind  conditions.  In  addition,  instrumental  errors 
appear  as  a fraction  of  the  full-scale  wind  setting  chosen, 
which  was  10  m sec-1  for  the  case  shown.  Recent  modi- 
fication to  the  instrumentation  have  reduced  this 
source  of  error. 

Figs.  11  and  12  compare  the  wind  speed  and  direction 
calculated  from  the  component  measurement  and  the 
same  quantities  measured  at  the  30-m  level  on  the 
tower.  All  recordings  are  10-sec  averages.  The  wind 
recordings  at  the  tower  are  shifted  40  s$c  in  time  to 
compensate  for  the  travel  time  from  the  triangle  to  the 
tower.  The  differences  in  the  wind  speed  magnitude 
are  expected  because  of  the  difference  in  height  above 
ground.  Using  the  logarithmic  wind  profile  to  calculate 
the  wind  speed  at  the  30  m level  from  the  wind  speed 
measured  with  the  optical  technique,  we  obtain  the 
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Fto.  11.  Comparison  of  wind  speed  measured  at  the  30  m level 
on  the  tower  (dotted  line)  with  the  wind  speed  calculated  (torn 
the  triangle  measurement  (solid  line).  The  dashed  line  is  the 
calculated  logarithmic  wind  protile  using  a roughness  parameter 
■(ill  4i  m,  and  a friction  velocity  *a  of  20  cm  sec 


dashed  line  curve  shown  in  Fig.  11.  A roughness  param- 
eter *0  and  a friction  velocity  m«  of  4 cm  and  20  cm  sec  \ 
respectively,  were  used.  These  values  are  regarded  as 
reasonable  for  the  type  of  ground  caver  and  wind  speed 
we  encountered. 

The  a>r relation  between  the  two  types  of  measurt 
ments  is  very  good.  The  larger  variance  of  the  wind 
measured  at  the  tower  is  cvpected  because  this  is  a 
point  measurement,  while  the  optical  wind  measure- 
ment is  an  average  both  in  time  and  space. 


Fto.  12.  Comparison  of  wind  direction  measured  at  the  30  m 
level  on  the  tower  (doited  line)  and  that  calculated  from  the 
laser  wind  measurement  v solid  line). 


5.  Conclusion 

It  has  been  shown  that  a wind  measuring  system 
using  three  se|>ar.ite  optical  systems,  each  consisting 
of  a 3 m\\  lie  Ne  laser  and  two  photodiode  receivers, 
forming  a equilateral  triangle,  is  successful  in  measuring 
the  divergence  over  a small  enclosed  area.  By  using  two 
'of  the  optical  systems,  we  can  get  a space-averaged 
horizontal  wind  vec  tor.  The  thermal  plumes  as  seen  by 
the  acoustic  sounder  and  the  convergence  compare  well. 
The  correlation  between  the  divergence  and  the  vertical 
wind  measured  on  the  tower  shows  that  the  entrainment 
into  plumes  is  promotional  to  the  vertical  velocity. 

During  jx-riocls  with  large  convective  activity  a 
certain  periixlicity  in  the  development  of  plumes  is 
found. 
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DISPERSION  AND  SPECTRA  OF  GRAVITY  WAVES  PROBABLY 
GENERATED  BY  A CONVECTIVE  STORM 


By 

A G Kjelaas,  E E Gossard,  J M Young  and  W R Moninger 


ABSTRACT 


In  this  paper  wc  describe  a case  stiuly  of  progressive  gravity  waves  in  the  troposphere.  We  deduce  their 
propagation  speed  anil  direction,  and  their  dispersive  characteristics,  from  surface  spectra  of  wind  and 
pressure.  We  compare  the  results  with  direct  measurements  of  the  waves  made  over  an  array  of  spaced 
pressure  tensors  using  cross  spectrum  techniques  and  an  algorithm  developed  by  one  of  us  (Young). 
The  agreement  is  very  satisfactory  and  demonstrates  the  feasibility  of  monitoring  atmospheric  wave 
conditions  in  almost  real  time  using  only  sensors  (high  quality)  in  use  at  standard  weather  observation 
stations. 

The  waves  in  the  case  we  describe  were  highly  dispersive  and  apparently  were  generated  by  a violent 
thunderstorm  some  100  km  west  of  the  site  of  observation. 


INTRODUCTION 

Atmospheric  gravity  waves  and  convective  activity  are  important  contributors  to  the 
atmospheric  pressure  spectrum  in  the  band  of  periods  from  about  1 minute  to  1 hour 
(Gossard,  I960).  The  form  of  the  spectrum  of  displacement  of  the  density  surfaces  (or 
isentropes),  and  the  dispersion  characteristics  of  propagating  disturbances  reveal  much 
about  the  mechanism  of  their  generation  (Gossard,  1974).  In  this  paper  we  measure  the 
spectrum  and  the  dispersion  of  gravity  waves  apparently  generated  by  a severe  convec- 
tive storm  on  8 August  1972  in  southeastern  Colorado.  The  wave  spectra  and  the 
dispersion  of  the  waves  we  will  discuss  are  dramatically  different  from  the  spectra  and 
dispersion  of  waves  analyzed  by  Gossard  and  Sweezy  (1974)  which  they  concluded  were 
generated  by  dynamic  instability  associated  with  shear  in  the  middle  troposphere. 

A method  devised  by  Gossard  and  Sweezy  (1974)  will  be  used  for  determining  the 
dispersive  characteristics  of  the  observed  waves,  and  we  will  compare  it  with  direct 
measurements  of  the  speed  and  direction  of  wave  components  across  an  array  of  spaced 
pressure  sensors.  For  the  latter  analysis  we  use  the  digital  techniques  and  computer 
algorithm  devised  by  one  of  us  (Young)  to  carry  out  a cross  spectral  analysis  between 
spaced  sensors  in  such  a manner  that  the  "beam"  of  the  array  is  effectively  steered,  or 
swept,  continuously  through  360  deg  of  azimuth. 

In  addition  to  the  usual  temporal  filtering  commonly  used  to  extract  mesofrequency 
information,  we  have  used  sensors  capable  of  spatially  filtering  the  atmospheric  spect- 
rum, and  in  this  paper  we  will  discuss  the  merits  of  such  spatial  filtering  for  sorting  out 
the  contribution  of  atmospheric  waves  for  analysis.  We  will  show  that  by  use  of  such  a 
spatial  averaging  of  the  wind  field,  the  large  mesoscale  features  in  the  wind  spectrum  can 
be  effectively  extracted  from  the  turbulence  ’’noise". 


2 ANALYTICAL  TECHNIQUES 


In  the  method  of  Gossard  ami  Sweezy  (1974),  the  spectra  of  surface  pressure  and  wind 
are  combined  to  obtain  the  "spectrum”  of  intrinsic  phase  velocity  of  waves  in  the 
mid- frequency  band  of  the  spectra  of  atmospheric  properties.  Within  the  constraints  of 
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the  linearized  theory  and  the  Boussinesq  approximation,  they  point  out  that  (c  g,  see 
Gossard  and  Hooke,  1974) 


- ?-jr0u  <>■*> 

In  these  equations 

= IP0(z)/^sl,/Ju>T?  and  P = Ip0(z)/Psr,/Jp 

where  u,  p,  and  r?  are,  respectively,  wave  perturbations  of  the  horizontal  wind,  pressure 
and  isentropic  surfaces;  g0(z.)  is  the  unperturbed  density  and  ps  is  the  density  at  some 
reference  level  such  as  the  earth's  surface.  Ihe  wave  frequency  relative  to  an  observer 
moving  with  the  mean  flow  is  co  (the  intrinsic  frequency)  and  the  wavenumber  k = 2tr/X 
where  X is  wavelength.  Then  the  displacement  spectrum  at  height  11  within  a surface 
layer  of  constant  Vaisala  Brunt  frequency  N and  constant  wind  is  related  to  the  pressure 
spectrum  at  the  surface  by 
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where  B is  sin  n|  H/n!  H or  sin  h !!/>]  II  depending  on  whether  Nt  /co  is  greater  than 
or  less  than  unity.  The  dispersion  equation  for  waves  of  infinitesimal  amplitude  relating 
n to  the  wave  frequency  co  and  wavenumber  k through  the  Vaisala-Brunt  frequency 
^=V/(g /#)  dO/dz  (where  0 is  the  potential  temperature)  is  essentially  (see  Gossard  and 
Hooke.  1974) 


n = i> 


(3) 


for  waves  whose  vertical  wavelength  is  small  compared  with  the  atmospheric  scale 
height. 


Also,  the  relation 
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provides  a measure  of  the  wave  dispersion  when  the  spectra  contain  an  important  wave 
contribution.  The  wind  component  U must  be  extracted  vcctorially  from  fluctuations  of 
the  total  wind  its  illustrated  in  figure  1.  I his  is  accomplished  by  choosing  a suitable 
numerical  low-pass  filter  to  apply  to  wind  speed  and  direction,  thus  separating  the  mean 
wind  from  the  total  wind,  lhe  same  low-pass  filter  is  applied  to  the  pressure  record. 
However,  our  pressure  sensing  system  had  its  own  high-pass  filter,  with  a time  constant 
of  50  sec,  whose  effect  had  to  be  removed  from  the  final  spectrum.  A low-pass 
numerical  filter  was  chosen  such  that  a perturbation  of  15  min  period  would  have  its 


Hgure  1 1 1'dor  diagram  Jor  extracting  the  wind  component  l’  along  the  ware  from 

fluctuations  of  the  total  wind 
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amplitude  reduced  to  one-third  the  true  value.  The  time  series  of  wind  and  pressure  were 
sampled  every  second  and  then  averaged  over  20  sec  intervals. 

Young’s  method  of  processing  data  from  spaced  sensors  is  to  azimuthally  steer  the  array 
t>eam  in  either  the  time  or  frequency  domains.  It  differs  from  other  algorithms  such  as 
that  of  Mach  and  Smart  (1972)  primarily  in  having  greater  flexibility  in  choice  of 
bandwidth  when  signal  velocity  is  approximately  independent  of  frequency  as  in  the 
case  of  acoustic  waves.  This  is  accomplished  by  processing  in  terms  of  the  wave 
’’slowness”  (defined  as  the  inverse  of  velocity)  rather  than  processing  in  frequency-wave- 
number space.  Varying  the  frequency  bandwidth  has  little  influence  on  the  average 
within  the  band  when  "slowness”  is  relatively  constant  over  the  band,  so  the  algorithm 
is  adaptive  in  choosing  an  optimum  trade-off  between  confidence  and  spectral  resolu- 
tion. When  applied  to  gravity  waves,  the  assumption  of  little  dispersion  is  no  longer 
necessarily  valid  as  we  shall  see  when  results  are  discussed. 


3 OBSKRVATIONAI.  TECHNIQUES 

The  pressure  sensors  used  in  this  study  were  essentially  like  that  described  by  Bedard 
(1971),  except  that  the  filter  amplifiers  producing  the  various  bandpasses  shown  in  his 
figure  3 were  not  employed.  Therefore  the  high  frequency  response  was  practically 
unlimited  (as  far  its  gravity  waves  are  concerned)  and  the  high-pass  characteristic  can  be 
described  by  the  single  time  constant  of  50  sec.  Each  presssure  sensor  had  30  m of  hose 
attached,  into  which  hypodermic  needles,  used  its  orifices,  were  inserted  at  intervals  of 
1.5  m.  I he  hose  serves  to  spatially  filter  out  small  scale  turbulent  inhomogeneities  in  the 
pressure  field.  The  total  pressure  array  consisted  of  a small  array  within  a larger  array  as 
shown  in  Figure  2.  The  small  array  was  designed  primarily  to  look  at  relatively  small 
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hgure  2 Layout  of  the  experimental  configuration 

The  solid  triangle  is  the  laser  triangle.  Points  marked  Ml,  M2,  MS  are  the  small 
mil  rnharograph  array,  while  M4,  M!>,  Mt>,  M7  denote  the  huger  array.  M7  is  20  m from 
thr  tower  base. 
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scale  features  within  the  boundary  layer,  and  the  large  array  was  designed  to  observe 
larger  scale  features  such  as  infrasound  or  long  waves  higher  in  the  atmosphere.  In  the 
present  study,  additional  spatial  filtering  of  the  pressure  field  was  achieved  by  averaging 
pressure  over  the  small  array  of  sensors  to  obtain  the  spectrum  of  pressure  for  use  in  Eqs 
(2)  and  (4). 

The  anemometers  used  in  the  present  study  were  bivanes  with  a high  frequency  response 
limited  to  approximately  one  Hz.  However,  our  primary  source  of  wind  data  was  a 
unique  system  designed  and  built  by  G Ochs  and  more  completely  described  by 
Lawrence  et  al  (1972)  and  Kjelaas  and  Ochs  (1974).  The  system  uses  a laser  path  to 
measure  the  average  wind  (somewhat  weighted  toward  the  center  of  the  path)  transverse 
to  the  path  by  measuring  the  velocity  of  the  scintillation  pattern  at  the  receiver  end. 
Clearly,  only  two  such  paths  are  needed  to  obtain  a value  of  the  total  horizontal  wind. 
However,  the  system  was  used  in  a triangular  configuration  so  that  horizontal  mass 
convergence  and  divergence  could  be  measured  (Kjelaas  and  Ochs,  1974).  This  is  also  a 
suitable  configuration  for  obtaining  wind  spatially  averaged  over  the  area  of  the  triangle. 
The  laser  triangle  is  shown  in  Figure  2. 

I he  results  of  the  analysis,  using  the  spatially  averaged  pressure  and  winds  as  described 
above,  were  compared  with  results  using  point  sensors.  The  advantage  of  using  the  laser 
winds  was  especially  evident,  and  values  of  U obtained  from  bivane  and  front  the  laser 
triangle  arc  shown  for  comparison  in  Figure  4.  Apparently  the  removal  of  small  scale 
turbulence  "noise"  by  the  spatial  filtering,  accomplished  by  laser  path  averaging  of  the 
wind  field,  significantly  improves  the  wave  signal  to  noise. 


4 EXPERIMENTAL  RESULTS 

The  case  chosen  for  analysis  was  an  event  which  occurred  in  the  early  morning  hours 
(0100— 0230  MI)T)  on  8 August  1972  at  Haswell,  Colorado.  Records  of  the  waves  of 
this  event  have  been  analyzed  using  an  acoustic  sounder  (mangle  and  the  results 
published  by  Kjelaas  et  al  (1974)  from  whose  paper  the  acoustic  sounder  records  in 
figure  3,  showing  most  of  the  time  interval  we  analyzed,  have  been  extracted.  It  is  the 
purpose  of  our  paper  to  analyze  the  dispersion  and  spectral  characteristics  using  the 
observational  techniques  described  above  and  try  to  deduce  the  mechanism  of  wave 
generation. 

The  pressure  record  and  the  unfiltered  and  low-pass  filtered  wind  speed  from  the  laser 
triangle  arc  shown  in  F'igure  4.  The  results  of  carrying  out  the  vector  subtraction 
indicated  in  Figure  1 arc  shown  in  Figures  5 and  6.  F'igure  5 shows  the  time  series  of  the 
component  U from  the  laser  winds  along  with  the  corresponding  quantity  from  the 
bivane  at  30  m elevation  on  a nearby  tower.  The  fluctuating  component  due  to  the 
waves  is  much  more  evident  in  the  laser  wind  measurements.  We  have  chosen  U as 
positive  for  clockwise  rotation  of  the  wind  vector.  Therefore,  U and  P are  positively 
correlated  when  the  waves  come  from  the  180  deg  sector  to  the  left  of  the  mean  wind, 
but  they  are  negatively  correlated  when  they  come  from  the  right-hand  180  deg  sector 
(as  is  the  case  for  the  principal  peak  in  the  directional  distribution  shown  in  Figure  6). 
Some  phase  shift  results  front  the  high-pass  filtering  of  the  pressure  fluctuations  by  the 
pressure  instrument.  Figure  6 shows  the  directional  distribution  of  the  waves  found  by 
the  method  of  Gossard  and  Sweezy  (1974).  The  arrows  superimposed  on  the  direction 
histogram  show  the  directions  obtained  using  Young’s  algorithm  and  the  spaced  array 
data.  The  arrow  from  270  deg  applies  to  the  band  from  21  min  to  80  min.  and  that  from 
225  deg  applies  to  the  entire  band  from  4 to  80  min. 

Figure  7 shows  the  dispersion  plot  of  phase  velocity  vs  frequency  by  the  method  of 
Gossard  and  Sweezy,  together  with  the  results  of  Young's  analysis  from  spaced  array 
data  shown  for  frequency  bands  indicated  by  the  bars.  The  spacing  of  elements  in  the 
array  produces  a high-frequency  cut-off  in  the  array  spectra.  The  excellent  agreement 
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Figure  3 Acoustic  sounder  record  of  part  of  the  event  analyzed  in  this  paper 

LT  is  Mountain  Daylight  Time.  To  convert  to  UT  add  (>  hours  to  L.T.  (From  Kjelaas  et  at 
(1974)). 


Figure  4 Pressure  recording  together  with 
the  low- pass  filtered  and  unfiltered 
laser  wind 

No  scale  is  included  for  pressure  since  the 
instrument's  filter  characteristic  has  not 
been  removed. 


Figure  3 Wind  component  V along 
the  ware  direction 
Upper  curve  is  obtained  using 
the  laser  triangle,  while  lower 
curve  is  from  the  bivane  at 
SO  m on  a nearby  tower.  Power 
density  of  frequencies  greater 
than  ~l?xl03  II/.  is  about  two 
orders  of  magnitude  down 
from  the  peak  in  the  power 
spectrum,  so  such  points  pro- 
bably have  no  significance. 
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is  a very  convincing  argument  for  the  validity  of  both  techniques.  At  a frequency  of 
2 * 10'3  Hz,  the  power  in  the  spectrum  (see  Figure  10)  has  diminished  by  two  orders 
of  magnitude  from  the  peak  power.  Thus  we  do  not  consider  the  velocities  in  Fi- 
gure 7 to  be  significant  at  frequencies  higher  than  about  2 x 10'3  Hz.  During  this 
wave  event,  the  average  surface  winds  were  very  light  (about  l ms’1 ) and  were  mostly 
between  north  and  east  in  opposition  to  the  wave  direction.  We  have  attempted  to 
obtain  C,  the  true  wave  velocity  relative  to  the  ground,  by  subtracting  1 ms'1  from 
oj/k.  This  actually  lessened  the  agreement  in  the  results  of  the  two  techniques. 


f(Hz) 


Figure  7 Phase  velocity  vs  frequency 

The  horizontal  bar*  are  the  phase  velocities  for  different  frequency  band*  found  by 
applying  Young’s  technique  to  the  large  microbarograph  array. 
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Figure  8 shows  the  plot  of  to  vs  k.  On  this  plot  the  slope  of  the  line  from  the  origin  to 
any  point  on  the  eurve  represents  the  intrinsic  phase  velocity,  and  the  slope  of  the 
tangent  to  the  curve  is  the  group  velocity  given  by  dco/dk. 


Figure  8 Intrinsic  frequency  as  a function  of  the  wavenumber 
t)nly  the  fundamental  mode  is  shown  for  model  no.  1. 

Figure  9a,b,c  shows  the  RAW1NSONDE  data  obtained  at  0457  MDT  (the  sounding 
nearest  in  time  to  the  wave  event).  Unfortunately,  the  sounding  does  not  well  represent 
the  lowest  layers  at  the  time  of  the  waves,  because  a high  surface  wind  sprang  up  at 
about  0400  MDT  apparently  due  to  outflow  from  the  cool  downdraft  of  a thunder- 
shower occurring  locally  at  about  that  time.  The  high  wind  at  the  surface  mixed  the  air 
within  the  radiation  inversion  and  elevated  it.  We  have  used  the  surface  wind  and 
temperature  data  at  the  time  of  the  wave  event  and  the  temperature  and  wind  data  at 
the  top  of  a nearby  151  tn  tower  (see  Figure  2)  to  attempt  to  correct  the  sounding  data 


Figure  ft  Radiosonde  at  0457  MPT 

I he  circles  shown  are  the  wind  and  temperature  data  from  the  tower  at  the  time  of  the 
wave  event.  The  various  atmospheric  models  discussed  are  shown  as  dashed  curves. 
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Figure  10  Spectrum  of  layer  displacement  at  a 
height  of  200  m 


back  to  the  time  of  the  event.  The 
values  arc  shown  as  circles,  and 
two  atmospheric  models  we  ana- 
lyze in  the  next  section  are  shown 
as  dashed  and  dash- dot  curves  in 
Figure  9c. 

Figure  10  shows  the  spectrum  of 
layer  displacement  at  a height  of 
200  m calculated  from  equation 
(2).  Both  the  spectrum  in  Figure 
10  and  the  dispersion  in  Figure  8 
differ  significantly  from  the  cor- 
responding quantities  found  by 
Gossard  and  Sweezy  from  the  data 
taken  at  San  Diego,  Calif,  under 
conditions  of  wind  shear  and 
instability  aloft.  The  spectrum  in 
Figure  10  implies  an  rms  ampli- 
tude of  38  m for  the  waves  at  a 
height  of  200  m. 
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5 COMPARISON  WITH  MODELS 


The  light  winds  aloft  shown  in  Figure  9 and  the  relatively  low  surface  and  tower  winds 
at  the  time  of  the  event  (compared  with  the  wave  velocity  of  about  10  ms'1  at  the 
spectrum  maximum)  imply  that  shear  is  not  the  generating  mechanism  for  these  waves. 
Furthermore,  the  waves  arc  propagating  against  the  background  winds.  Therefore  it 
seems  reasonable  to  compare  the  predictions  of  a shearless  model  with  the  observed 
dispersion. 


We  adopted  a shearless  model  consisting  of  three  layers  rather  than  more  realistic  models 
with  many  layers  because  a three-layer  model  can  be  studied  analytically  whereas 
models  with  four  or  more  layers  are  more  effectively  studied  using  numerical  methods 
on  high  speed  computers.  We  used  a normal  mode  approach  leading  to  an  eigenvalue 
equation  whose  derivation  is  given  by  Gossard  et  al  (1970).  They  find  that 


n2  cot  2n2  All  = 


n2a-n,  y3  cot  n,H 
n,  cot  n(H  + y3 


(5) 


where  the  subscripts  apply  to  the  three  layers  of  a model  shown  schematically  in 
Fipire  12.  Equation  (5)  was  applied  to  the  two  models  shown  by  the  dashed  and 
dash-dot  curves  in  Figure  9c  (see  caption  of  Figure  12).  In  model  no.  1 the  thickness 
2AH  of  the  middet  layer  is  zero  so  it  is  a two-layer  model  for  which  equation  (5) 
becomes  (see  Gossard  and  Hooke,  1974) 
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small  for  the  wavenumber  range 
we  consider  here.  The  waves  stu- 
died by  Gossard  and  Sweezy 
(1974)  displayed  the  non-disper- 
sive  character  to  be  expected  of 
wave  generation  by  shear.  The 
highly  dispersive  nature  of  the 
waves  we  have  studied  in  this 
paper  implies  a different  mecha- 
nism of  generation.  In  addition. 
Figure  9 shows  that  the  small 
shear  that  might  exist  during  the 
period  of  the  event  would  lead  to 
waves  propagating  from  the  north- 
easterly direction.  They  were 
actually  propagating  from  the 
west  through  south-west  sector  (as 
seen  in  Figure  6)  with  a group 
velocity  of  about  6 ms'1  at  fre- 
quencies corresponding  with  the 
maximum  in  energy  of  the  spec- 
trum shown  in  Figure  10.  If  we  look  100  km  to  the  west  of  Haswell  on  the  radar 
weather  map  shown  in  Figure  1 1 taken  at  2135  MST  (2235  N1DT)  obtained  by  the  NWS 
using  the  WSR  57  radar  at  Liman,  Colorado,  we  find  a violent  thunderstorm  near 
Pueblo,  Colorado,  extending  above  39  000  ft  (11.9  km).  By  2337  MDT  the  tops 
exceeded  41  000  ft  (12.5  km).  It  seems  very  likely  that  the  waves  discussed  in  this  paper 
were  generated  by  the  storm.  The  precise  mechanism  is  uncertain.  It  seems  possible  that 
the  storm  may  act  like  a slow  explosion  in  the  atmosphere  effectively  generating  the  low 
frequency  gravity  waves  in  our  example.  In  that  case  the  analysis  of  Ramm  and  Warren 
(1963)  of  gravity  wave  generation  by  an  impulsive  source  may  be  relevant.  Note  that 
Figure  10  shows  most  of  the  energy  to  be  centered  at  a wave  period  of  about  14  min. 
On  the  other  hand,  it  is  very  possible  that  the  cool  gust  front  associated  with  the  strong 
downdraft  of  the  storm  generated  the  waves  in  the  manner  of  a travelling  disturbance. 


6 CONCLUSION 

We  conclude  that  the  methods  described  in  this  paper  and  that  of  Gossard  and  Sweezy 
(1974)  provide  a useful  technique  for  observing  the  properties  of  waves  in  the  atmos- 
phere and  of  diagnosting  the  probable  mechanism  of  wave  generation.  The  method 
requires  only  sin^e  station  measurements  of  pressure  and  wind.  However,  comparison  in 
Figure  5 of  the  wind  spatially  averaged  over  the  laser  triangle  with  measurements 
obtained  with  a point  sensor  emphasizes  the  great  value  of  spatial  filtering  in  the 
extraction  of  large  (or  meso)  scale  features  from  the  turbulence  dominated  total  wind 
field  spectrum. 
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Figure  12  The  three-layer  model 

Model  no.  1 (tlash  dot  curve  in  Figure  8c)  is 
a two-layer  model,  i e.  Ail  = 0,  in  which 
Ni  - 0.0S54,  N,  = 0.01 1 1 and  II  = 200  m. 
In  model  no  2 (dashed  curve  in  Figure  8c). 
N | - 0.050,  N2  « 0.0202.  N,  = 0.01 11, 
H = 100  in.  2 All  - 600  m. 
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WAVES  OBSERVED  IN  THE  PLANETARY  BOUNDARY  LAYER 
USING  AN  ARRAY  OF  ACOUSTIC  SOUNDERS 

A G Kjelaas 

D W Beran  , W II  Hooke  , and  B R Bean 


I ABSTRACT 

Three  acoustic  sounders,  with  antennas  having  vertically-pointed  beams  located  at  the  vertices  of  a 
triangle  about  300  m on  a side,  have  been  used  successfully  to  calculate  horizontal  phase  velocities 
, of  gravity  waves  in  the  lowest  1000  m of  the  atmosphere.  The  wave  parameters  obtained  were 

compared  with  those  obtained  by  an  array  of  microbarographs.  The  results  from  the  two  indepen- 
dent techniques  agreed  well  in  all  cases  where  the  waves  were  detected  by  both  methods.  However, 

Ithe  acoustic-sounder  array  was  also  able  to  detect  gravity  waves  propagating  in  an  inversion  over- 
lying the  convective  boundary  layer  that  were  undetected  by  the  microbarograph  array  because  of 
wave  attenuation  in  the  convective  region. 

Wave-associated  vertical-velocity  fluctuations  inside  the  undulating  structure  were  measured  using 
Doppler  techniques  to  an  accuracy  ot  about  f 0.1  m sec'1 . Vertical-velocity  fluctuations  of  about 
-0.5  m sec  1 were  measured  inside  one  wave  with  displacement  amplitude  of  120  m and  wavelength 
of  5.5  km.  From  wave-associated  vertical-velocity  fluctuations,  displacement  amplitude,  and  wave 
period,  we  could  calculate  the  wind-speed  component  along  the  wave  propagation  direction.  This 
calculated  wind-speed  component  agreed  reasonably  well  with  the  value  derived  independently  from 
radiosonde  measurements. 


1.  Introduction 

During  the  summer  of  1972,  the  staff  of  NOAA’s  Wave  Propagation  Laboratory 
(WPL)  mounted  an  experimental  program  at  the  Haswell  field  site  in  Colorado  lasting 
several  weeks.  During  the  observing  period,  this  area  was  characterized  by  a typical 
High  Plains  meteorological  regime:  strong  radiative  cooling  at  night  and  strong  heating 
and  convective  activity  during  the  day.  The  program  had  the  dual  objectives  of  (1) 
checking  the  operation  of  a variety  of  remote  sensing  systems  and  calibrating  them 
against  data  gathered  from  a 150-m  instrumented  meteorological  tower,  and  (2) 
studying  the  Haswell  atmospheric  boundary  layer. 

Included  among  the  remote-sensing  systems  and  meteorological  instruments  were 
three  acoustic  echo  sounders  and  an  array  of  microbarographs.  The  acoustic  echo 
sounders  operate  like  sonars;  they  transmit  acoustic  pulses  and  detect  the  echoes 
resulting  from  boundary-layer  turbulent  wind  and  temperature  fluctuations  (Tatarski, 
1961;  Monin,  1962).  In  this  way,  these  sounders  can  measure  boundary-layer  turbu- 
lence intensity  (eg  Little,  1969;  Bean,  1971;  Emmanuel  et  al,  1972;  Beran  et  al, 
1973;  McAllister,  1968)  as  a function  of  height  and  time,  as  well  as  the  three  compo- 
nents of  the  boundary-layer  wind,  which  can  be  determined  from  the  Doppler  shift  of 
the  acoustic  echoes  when  two  of  the  three  antenna  beams  are  directed  at  angles  off 
the  vertical  (Beran  et  al,  1971,  1973;  Beran  and  Clifford,  1972).  The  microbarographs 
supported  the  remote-sensing  studies  of  boundary-layer  gravity  waves,  by  measuring 
wave-associated  pressure  fluctuations  at  the  earth’s  surface.  From  these  records  we 
determined  wave  parameters  such  as  wave  period,  horizontal  wavelength,  phase  velo- 
city and  azimuth  of  propagation  (for  illustrations  of  previous  studies  see  , e g,  Hooke 
etal,  1972,  Beran  et  al,  1973).  

During  most  of  the  observing  period,  the  acoustic  echo  sounder  data  were  gathered 
primarily  to  evaluate  the  relative  merits  of  various  antenna  configurations  in  Doppler 
wind  measurements.  However,  we  were  also  interested  to  see  whether  we  could  ope- 
rate the  sound  is  in  a mode  analogous  to  that  of  the  microbarograph  array  in  order 
to  study  boundary-layer  wave  motions.  To  test  this  possibility,  we  operated  the 
sounders  for  short  intervals  several  hundred  meters  apart  with  all  beams  directed 
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vertically.  Then  to  determine  horizontal  wavelengths,  wave  phase  speeds,  and  propaga- 
tion directions,  we  used  the  cross-correlation  between  corresponding  features  of  the 
acoustic  echo  returns  recorded  by  the  three  sounders  (i  e,  the  cross-correlation  be- 
tween the  height- time  echo  intensity  patterns),  and  the  time  lags  required  to  maxi- 
mize this  cross  correlation.  These  values  were  compared  with  those  obtained  from  the 
microbarograph  cross-correlation  analysis  as  well  as  with  those  obtained  by  using  the 
so-called  impedance  relation  (Gossard  and  Munk,  HIM),  which  gives  the  wave  phase 
speed  and  direction  from  an  analysis  of  the  wave-associated  pressure  and  velocity 
fluctuations  measured  at  a single  point. 

In  this  paper,  we  summarize  the  results  of  using  an  array  of  three  vertically-pointing 
acoustic  echo  sounders  in  calculating  phase  velocity,  period,  and,  hence  also  horizon- 
tal wavelength  of  gravity  waves.  The  results  illustrate  the  advantages  and  limitations 
of  this  application  of  acoustic  echo  sounding.  In  addition,  we  show  that  from  sounder 
measurements  of  (l)wave  period,  (2)  wave- associated  vertical-velocity  perturbations 
(as  revealed  by  Doppler  shifts  of  the  acoustic  echo  returns),  (3)  vertical  parcel- 
displacement  perturbations  (as  revealed  by  fluctuations  in  height  of  the  strata  of 
strong  echo  returns),  and  (4)  horizontal  wavelength  (calculated  from  measurements  of 
wave  period  and  phase  velocity),  we  can  determine  the  mean  wind  speed  component 
in  the  direction  of  wave  propagation.  We  can,  therefore,  check  the  data  still  further 
for  consistency  by  comparing  the  value  of  the  mean  wind  speed  calculated  in  this 
manner  with  that  measured  directly  from  a radiosonde  and  on  the  tower.  Thus  we 
have  demonstrated  in  principle  that  acoustic  echo  sounders  can  use  boundary-layer 
gravity  waves  as  tracers  in  wind  profiling. 


2.  The  Kxperiment 


The  experimental  configuration  is  shown  in  Fig  1.  The  points  marked  A,  B,  and  C arc 
the  positions  of  the  three  acoustic  echo  sounders,  which  form  a right  triangle  with 
legs  roughly  300  m in  length.  1'he  points  marked  Ml,  M2,  and  M3  denote  the  posi- 
tions of  the  three  microbarographs  and  T the  location  of  the  150  m meteorological 
tower.  This  size  and  configuration  was  adequate  for  the  gravity-wave  studies  here; 
however,  to  study  larger-scale  waves  with  higher  phase  velocities,  we  would  have 
needed  a larger  triangle. 
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PEAK  POWER 

10  W(  ACOUSTIC 

PULSE  WI0TH 

100  m sac 

PULSE  REPETITION 
FREQUENCY 

0.126  sac  1 

MAXIMUM  RANGE 

1380  m 

CARRIER  FREQUENCY 

VARIABLE  FROM 
1760-2780  Hi 

ANTENNA  DIAMETER 

162  cm 

ANTENNA  BEAMWI0TH 

±6°  TO  3dB  POINTS 

BEAM  DIRECTION 

VERTICAL 

RECEIVER  BANDWI0TH 

100  Hi 

Table  1 Acoustic  Echo  Sounder  Parameters 


I'hc  three  sounders.  A,  B,  and  C, 
were  operated  at  center  frequencies 
of  2750  Hz,  1 750  Hz,  and  2250  Hz. 
respectively,  in  order  to  avoid  cross 
talk.  The  other  sounder  parameters 
are  given  in  Table  i.  Both  the  inten- 
sity of  the  backscattercd  energy  and 
the  Doppler  frequency  shift  for  each 
sounder  were  recorded  on  analog- 
tapes  and  digitized  at  a sampling  rate 
of  62.5  sec'* . In  addition,  facsimile 
recordings  were  made  of  the  acoustic- 
echo  returns  as  functions  of  boun- 
dary-layer height  and  time.  The 
tower  was  instrumented  to  measure 
wind  and  temperature  at  30  m, 
60  m,  90  m,  120  m,  and  150  m;  and 
to  measure  wave-associated  pressure 
fluctuations  at  90  m and  150  m. 
Radiosonde  wind  and  temperature 
profiles  were  also  obtained. 


3.  Theory  of  Wind  Speed  Determination  from  Gravity-Wave  Observations 


for  any  fluid  motion  the  vertical  parcel  velocity  w is  instantaneously  related  to  the 
vertical  parcel  displacement  f by 


w 


= W 

Dt 


(1) 


where  D/Dt  denotes  the  time  derivative  billowing  the  parcel  motion.  That  is. 


D_  m 

Dt  dt 


+ u- 


V 


(2) 


so  that  for  a sinusoidal  wave  of  frequency  u>  (as  measured  by  an  observer  fixed  with 
respect  to  the  ground)  and  wave  vector  k,  and  of  perturbation  amplitude, 


(3) 


lhe  frequency  SI  is  the  so-called  "intrinsic”  wave  frequency,  i c,  the  wave  frequency 
that  would  be  measured  by  an  observer  moving,  relative  to  the  ground,  with  the 
background  wind  velocity  UQ.  In  general  UG=  UQ(z)  and  SI  = S2( z)  where  z is  the 
height.  Thus,  at  any  height  z,  we  have  the  relation 


w*(z)  = S2(z)  f*(z) 


H) 


where  the  asterisks  denote  amplitudes  of  the  wave-associated  fluctuations.  We  have 
then  that 


Uok(z)  = « 


w*(z) 

kr(*) 


r(1-2p> 


(-r>) 


This  value  of  U0^  can  then  be  compared  with  the  one  obtained  from  radiosonde. 


(> 


4.  Analysis 


HASWKIX.  C Ol.OKAllO 
August  8.  1972 


Three  periods  - 0100-01 40LT  Augusts,  1972  (Case  I);  021 5-0235LT,  July  23, 
1972  (Case  II);  and  0820-08501  ,T,  July  23,  1972  (Case  III)  - were  selected  for 

detailed  analysis.  The  facsimile 
records  from  each  of  the  three 
sounders  for  each  of  the  three 
cases  are  shown  in  Figs  2,  3, 
and  4.  To  determine  the  hori- 
zontal wave  phase  speed  and 
direction  of  propagation  of  the 
waves  revealed  by  the  undula- 
ting strata  of  strong  echo  re- 
turns shown  in  Figs  2,  3,  and  4, 
a simple  time-lagged  cross-corre- 
lation analysis  between  the  in- 
tensity patterns  of  the  strong 
echo  return  above  the  nocturnal 
inversion  for  the  three  sounders 
was  performed.  In  general,  the 
peak  in  the  cross  correlations 
does  not  occur  at  zero  time  lag, 
but  is  shifted  in  time.  These 
shifts  correspond  to  the  travel 
time  of  a wave  from  one  soun- 
der to  another.  From  this  travel 
time  both  phase  speed  and  di- 
rection of  propagation  can  be 
derived  (Briggs  et  id,  1950). 


0200 


0100 


0000 


The  accuracy  in  the  derived 
phase  velocity  depends  on  the 
accuracy  with  which  the  diffe- 
rent travel  times  can  be  estima- 
ted. Since  the  pulse  repetition 
frequency,  which  can  be  inter- 
preted as  the  sampling  frequen- 
cy, is  0.125  Hz  in  our  case,  the 
position  in  time  of  the  peak  of  the  cross  correlation  cannot  be  estimated  more 
accurately  than  ±4  sec.  1'his  should  give  an  accuracy  in  the  calculated  phase  speed  of 
the  order  of  ±1  m sec'1  and  approximately  ±10°  on  direction  for  a wave  travelling 
with  a speed  of  10  msec"1  over  the  array.  An  increase  in  phase  speed  will  increase 
the  estimated  errors.  For  Case  I,  Doppler  data  were  used  to  estimate  the  wave- 
associated  vertical-velocity  fluctuations.  For  this  purpose  the  echo  returns  from  five 
vertical  height  intervals  totalling  about  1 1 m were  time-averaged  over  10  consecutive 
pulses  (80  sec.).  The  vertical  velocity  resolution,  assuming  that  the  signal  must  be 
observed  for  0.0588  sec.  and  that  one-quarter  of  the  period  is  needed  to  track  the 
Doppler,  was  found  to  be  ±0.1  m sec'1. 


Local  Time 


Figure  2 Facsimile  records  from  each  of  the 
sounders  for  0000  02001.1'  August  8, 
1972:  Case  1. 


Case  1:  0100  01 401. T August  8,  1972.  In  Case  I (Fig.  2)  we  have  an  apparently  stable 
undulating  stratum  of  strong  echo  returns  well  separated  from  the  nocturnal  surface 
boundary-layer  inversion  below.  The  wave  motion  appears  to  be  nearly  sinusoidal,  and 
the  wave  parameters  for  this  time  interval  calculated  from  both  the  acoustic-echo- 
sounder  array  <lata  and  the  miernbarograph  array  data  arc  shown  in  Table  2.  Table  2 
also  gives  the  horizontal  wave  phase  speed  u>/k  and  propagation  direction  computed 
from  tower  data  on  the  amplitude  of  the  wave-associated  pressure  fluctuations  p’  and 
horizontal  velocity  fluctuations  u\  (where  prime  denotes  the  wave-associated  depar- 


WAVE  PARAMETERS 


ACOUSTIC 

SOUNDER 


MICROBAROGRAPH 

ARRAY 


"IMPEOANCE 

RELATION 


\ 


PERIOD 

716  sec 

717  sec 

HORIZONTAL  PHASE  SPEED 

7.6  nt  sec  1 

8.0  m sec  1 

7.7  m sec1 

HORIZONTAL  WAVELENGTH 

6.4  km 

6.7  km 

DIRECTION  OF  ARRIVAL 

245 

O 

o 

c* 

N 

AVERAGE  HEIGHT  OF  THE 

WAVE 

400  m 

AMPLITUDE  OF  THE  WAVE 

120  m 

table  2 hViiv  Parameters  for  tht  t’ rents  of  Case  / 
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ture  from  the  background  value 
k),  using  the  "impedance  rela- 
tion” (Gossard  and  Munk, 
1954) 


k Po(z)ll’(z) 


1 ok(Z) 


(6) 


No  significant  differences  ap- 
pear between  the  results  obtain- 
ed from  the  microbarograph  ar- 
ras and  the  acoustic  echo  soun- 
dei  array  as  far  as  wave  period 
and  horizont.il  phase  speed  are 
concerned;  both  give  a period 
of  about  715  seconds  and  a 
phase  speed  of  about  8 m sec’1 . 
The  difference  in  the  observed 
direction  of  arrival  is  larger, 
about  25°;  however,  because 
the  wavelength  is  an  order  of 
magnitude  greater  than  the 
dimensions  of  the  array,  the 
direction  determined  is  subject 
to  uncertainties  larger  than 
±10°. 


The  Doppler-derived  vertical 
velocities  for  Case  l are  shown 
in  f ig.  5.  The  lined  portions  in- 
dicate height  time  intervals  of 
positive  (i  e,  upward)  vertical 
velocity,  while  the  dotted  por- 
tions correspond  to  areas  with  negative  (i  e,  downward)  velocity.  The  heavy  line  is  the 
zero  vertical  velocity  houndarv,  the  *0.1  and  ♦-0.9  msec  1 contours  are  also  shown. 
For  clarity,  the  Isoundaries  of  the  stratum  of  strong  echo  return  (dotted  lines)  arc 
also  shown. 


Figure  d Facsimile  records  from  cat  h of  the 
sounders  Jo r 0200  0.1001  I July  2.1. 
I <>72:  Case  II 


I'he  Doppler  derived  vertical  velocity  for  the  descending  portion  of  the  wave  between 
0109  .uul  0I15I.T  was  found  to  average  0.25  msec1  with  a maximum  value  of 
0.55  msec  1 . between  0128  and  01291.1,  the  Doppler  method  gave  about 
0.85  m sec  1 with  a maximum  of  0.48  msec  1 . Ilus  agreement  is  very  good.  The 
low  sign  alt  o-no  i se  ratio  precluded  our  deriving  the  vertical  velocity  for  that  part  ol 
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Figure  4 Facsimile  records  from  each  of  the 
sounders  for  0800- 0900LT  July  23. 
1972:  CaselU. 


the  atmosphere  between  the 
stratum  of  strong  echo  returns 
and  the  nocturnal  inversion  be- 
low; therefore,  we  were  unable 
to  determine  whether  this  part 
of  the  atmosphere  was  oscilla- 
ting in  phase  with  the  wave  mo- 
tion above.  However,  by  exami- 
ning the  vertical  velocity  fluctu- 
ations at  the  top  of  the  noctur- 
nal inversion,  as  shown  at  the 
bottom  of  Fig  5,  we  see  that 
periods  with  significant  upward 
vertical  velocities  correspond  to 
the  periods  of  upward  vertical 
velocity  inside  the  stratum  of 
strong  echo  returns.  The  facsi- 
mile record  in  Fig  2 also  gives 
evidence  that  the  height  of  the 
nocturnal  inversion  is  oscillating 
in  phase  with  the  wave  motion. 

Using  the  wave  parameters 
shown  in  Table  2 together  with 
Eq  (5),  we  obtained  a wind 
component  of  6 m sec'1  along 
the  propagation  direction  of  the 
wave.  The  radiosonde  gave  a 
component  of  about  5.5  m 
sec'1.  Unfortunately,  the  radio- 
sonde we  used  was  released 


about  two  hours  before  the  wave  event.  During  the  time  between  the  radiosonde  and 
the  wave  event,  no  significant  changes  in  wind  speed  or  direction  were  observed  on 
the  nearby  tower.  This  may  indicate  that  the  radiosonde  data  were  representative  of 
values  prevailing  at  the  lime  of  the  event,  at  least  in  the  first  few  hundred  meters. 
Taking  into  account  the  difference  in  time  between  the  radiosonde  and  the  wave 
event,  the  agreement  between  the  wind  component  measured  by  the  radiosonde  and 
the  one  calculated  from  acoustic  sounder  data  is  quite  good.  This  appears  to  confirm 
the  potential  of  acoustic  echo  sounder  observations  of  boundary-layer  gravity  waves 
for  measuring  boundary-layer  wind  speeds. 


Case  11:  0200  02301. T July  23.  1972.  In  Case  11  (F'ig  3)  the  wave  is  clearly  within  the 
nocturnal  inversion.  The  wave  parameters  calculated  both  from  the  acoustic  sounder 
and  the  microbarograph  array  arc  shown  in  Table  3. 
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Figure  5 Doppler-derived  vertical  velocities  for 
the  case  shown  in  Figure  2. 

The  darker  line  shaded  portions  indicate  upward  vertical 
velocity,  while  lighter  dot  shaded  portions  correspond  to 
areas  with  downward  vertical  velocity.  T he  heavy  line  is  the 
zero  vertical  velocity  contour;  the  ±0.1  and  i 0.5  m sec'1 
contours  arc  also  shown. 

The  dotted  lines  are  the  boundaries  of  the  stratum  of  strong 
echo  return. 


WAVE  PARAMETERS 

SOUNDER 

PERIOD 

20S  mc 

HORIZONTAL  PHASE  SPEED 

5.2  m $*c  ' 

HORIZONTAL  WAVELENGTH 

1.1  km 

DIRECTION  OF  ARRIVAL 

271" 

AMPLITUDE  OF  THE  WAVE 

~ 50  m 

Table  4 Wave  Parameters  for  the  Events  of 

Case  III 


Thus  we  find  that  the  wave  amplitude  at  ground 
it  is  in  the  overlaying  inversion;  this  causes  it  to 
the  large  pressure  fluctuations  associated  with  the 


As  we  see,  an  excellent  agreement 
exists  between  the  two  systems.  In 
this  case  the  wavelength  is  much 
shorter,  only  about  3.5  km,  which 
gives  us  better  accuracy.  In  addi- 
tion, even  at  only  250  m above 
ground  the  wave  has  a relatively 
large  amplitude,  giving  us  a good 
signal-to-noise  ratio  in  the  pressure 
fluctuations  observed  on  the 
ground. 

Case  III:  0820  08451  T July  23, 
1972.  Case  III  (Fig  4)  is  a com- 
mon morning  situation  in  which 
the  nocturnal  inversion  is  elimina- 
ted when  a convective  layer  deve- 
lops with  the  onset  of  solar  heat- 
ing. The  rising  plumes  within  this 
convective  layer  appear  to  induce 
a wave  structure  at  its  top.  There 
is  some  evidence  of  vertical  wave 
phase  variation,  indicating  a wave- 
associated  vertical  energy  and 
momentum  transport  within  the 
overlaying  inversion.  These  wave 
parameters  are  shown  in  Table  4. 
For  this  event,  we  were  only  able 
to  calculate  the  wave  parameters 
from  the  acoustic  sounder  array; 
the  wave  was  not  detected  by  the 
microbarograph  array.  The  reason 
for  this  can  be  seen  from  the  gra- 
vity-wave dispersion  equation 
which,  in  the  Boussincsq  and  WKB 
approximations  can  be  written  as 
(c  g,  Hooke  et  al,  1973) 


nJ  = k’(g-l) 


where  we  have  assumed  that  the 
wave-associated  velocity  and  pres- 
sure perturbations  vary  as 

exp  i(Slt  - kx  - nz.) 

and  where  N*  = gd(ln  0)/dz  is  the 
square  of  the  Brunt  V'iiisala  fre- 
quency N.  Since  for  the  most  part 
the  temperature  structure  within 
the  convective  planetary  boundary 
layer  is  adiabatic,  NiO,  which 
implies  from  Kq  (7)  that  n = ± ik. 
level  is  only  a few  per  cent  of  what 
be  undetectable  in  the  presence  of 
thermal  plumes. 
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5.  Conclusion 

Wc  have  shown  that  three  vertically-pointing  acoustic  sounders  arranged  in  a triangle 
about  300  m on  a side  can  be  successfully  used  to  calculate  horizontal  wavelengths, 
phase  velocities,  and  direction  of  arrival  of  gravity  waves  in  the  lower  1000  m of  the 
atmosphere.  The  gravity  waves  analyzed  had  horizontal  wavelengths  ranging  from 
about  1 km  to  almost  6 km  with  horizontal  phase  velocities  from  5 to  10  msec'1. 
Comparing  these  results  with  a collocated  array  of  microbarographs  gave  a good 
agreement. 

Using  estimates  of  the  wave-associated  vertical-velocity  fluctuations  and  vertical 
displacements  together  with  the  wave  period  and  horizontal  wavelength,  we  were  able 
to  determine  the  mean  wind  speed  component  along  the  wave  propagation  direction. 
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ABSTRACT 

The  inthiencc  of  atmospheric  stability  on  the  vertical  spectra  of  retractive  index  and  air  velocity  is 
st  u diet! 

Spectra  of  refractive  index  and  air  velocity  are  deduced  from  (team  swinging  experiments  in  the  area 
txlwren  Norway  and  Denmark  where  a height  interval  from  some  400  m to  some  1*  km  is  involved  These 
spectra  are  compared  with  the  results  of  data  obtained  from  radiosondes. 
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1.  Introduction 

During  the  last  decade  a number  of  e\|>erinu'nts  have 
been  performed  from  which  definite  conclusions  have 
been  drawn  regarding  the  shape  of  the  refractive  index 
fluctuation  spectrum.  Manx-  of  these  results  are  in  verx 
poor  agreement  with  the  simple  turbulence  theory  based 
on  the  Kolmogorov  inertial  subrange  k y>  law.  Al- 
though it  is  probably  fair  to  say  that  some  of  the  earlier 
experiments  had  too  few  degrees  of  freedom  to  yield 
experimental  results  having  a significant  confidence 
level,  others  do  haxt*  a confidence  level  sufficient  to 
justify  a more  detailed  examination  of  tfic  findings.  It 
is  the  purpose  of  this  contribution  to  compare  the  re 
suits  from  simple  straightforward  radio  experiments 
with  results  from  conventional  direct  atmospheric 
soundings  in  the  hope  of  finding  a common  denominator. 

Specifically,  we  deduce  the  shape  of  the  refractive 
index  wavenumber  and  air  velocity  spectra  from  radio 
experiments  and  we  compare  these  with  data  obtained 
by  the  use  of  conventional  radiosondes. 

2.  The  radio  experiment 

The  radio  data  were  obtained  in  a bistatic  beam 
swinging  scatter  experiment  making  use  of  one  powerful 
transmitter  and  two  sensitive  receivers  located  far  he 
yond  the  horizon  relative  to  the  transmitter.  The  ex 
periment  was  performed  in  the  geographical  area  he 
tween  southern  Norway  (Kristiansand)  and  Denmark 
(Spodsbjerg  near  Copenhagen)  probing  the  atmosphere 
in  the  height  range  from  >■  5(K>  tn  to  ^ 11  kin. 

Ream-swinging  experiments  are  based  on  two  theo- 
retically obtained  equations,  the  first  of  which  is 

/>R0C  (1) 

where  /’«  is  the  received  scattered  power,  r(.r,y,j)  the 
homogeneity  function  relating  the  degree  to  which  the 


scattering  ability  (spectral  density)  varies  with  posi- 
tion (.v.v.sl  in  space,  and  <f>(  K)  the  refractive  index  spec- 
trum. The  vector  K is  measured  in  the  plane  in  which 
the  scattering  takes  place  such  that  K bisects  the  angle 
between  the  incident  and  the  scattered  wave  having 
wavenumbers  k,  and  k„,  respectively.  Accordingly,  K 
is  given  by 

K = k|—  k,  and  | K | = 2k  sin(0/2), 

where  t>  is  the  scattering  angle,  k = 2n  X,  and  X the  radio 
wavelength.  For  details  see  ('.jessing  (t‘)62). 

Whereas  the  first  basic  equation  relates  the  time- 
averaged  received  power  and  the  spatial  distribution  of 
refractive  index,  the  second  relates  the  frequency  of 
the  scattered  radio  waves  and  the  atmospheric  pattern 
of  motion.  Again  referring  to  ('.jessing  (1^021  for  details, 
we  have 

> 

(W  [I  sitvcf-MI'O)  sin(f>  2!]-  (21 

X 

where  SF  is  the  width  of  the  received  frequency  spec- 
trum (speed  of  fading),  l’  the  mean  cross  path  compo- 
nent of  the  wind  velocity,  i)  the  radio  beam  width,  51' 
the  fluctuating  component  in  the  direction  y (y  is 
essentially  the  direction  of  K relative  to  the  vertical 
direction),  and  0 the  scattering  angle. 

I he  objective  of  the  current  experimental  investiga- 
tion is  to  determine  <f>(K'  for  vertical  K from  (11  and 
to  determine  the  vertical  component  of  the  fluctuating 
velocity  51',  from  (2V  This  is  achieved  by  two  beam- 
swinging  experiments  (see  (1  jessing  and  Borresen, 
l%7). 

In  Experiment  A the  transmitter  and  receiver  beams 
are  elevated  sy  nchronously  in  the  ,v,  s great  circle  plane, 
thus  keeping  the  direction  of  K constant  and  vertical, 
whereas  K varies  since  t>  varies  (r  is  the  vertical 
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coordinate  whereas  a is  measured  along  the  path).  The 
received  |>ower  and  the  width  of  the  Doppler  spectrum 
(fading  spectrum)  are  measured,  From  this  experiment 
we  obtain  0( | K i ) front  (1)  for  K vertical  with  the  effect 
of  the  homogeneity  function  c(s)  superimposed  („•  being 
the  vertical  coordinate),  furthermore,  from  (2),  we 
obtain  the  mean  cross  path  wind  velocity  I and  the 
fluctuating  component  of  wind  velocity  51  along  the 
vertical  direction  since  K is  vertical. 

In  Experiment  B the  beams  are  also  swung  syn 
chronously  in  the  .v,s  plane,  but  in  this  case  in  such  a 
manner  that  the  elevation  angle  of  one  of  the  antennas 
increases  by  the  same  amount  that  the  elevation  angle 
of  the  other  decreases.  Hence,  the  scattering  angle,  and 
thus  | Kj  are  constant.  Note  that  the  direction  of  K is 
changed  through  the  small  angle  (some  3°' 

only  , whereas  the  position  of  the  scattering  volume  i-- 
lowered  and  shifter!  toward  the  receiver  (when  the 
elevation  angle  of  the  receiver  beam  is  increasing),  for 
all  practical  purposes,  therefore,  dir  k and  | K can  be 
considered  constant  while  x and  . are  varied.  It  is  well 
known  that  the  mean  refractive  index  varies  with 
height,  while  the  properties  are  generally  more  constant 
in  a horizontal  plane.  If  the  variation  in  e [Eq.  (1)]  is 
negligible  when  .r  varies  as  in  Experiment  B,  the  results 
of  this  experiment  can  be  used  to  compensate  for  the 
effect  of  atmospheric  nonhomogeneity  on  Experiment 
A.  By  performing  Experiment  B in  such  a way  that  the 
scattering  volume  is  shifted  to  either  side  of  the  mid- 
path position,  i.e.,  x is  varied  from  4-A.r  to  A.v,  while 
5 is  varied  twice  between  Ci  to  one  obtains  a measure 
of  the  influence  of  x on  >■  from  symmetry  considerations. 

Thus,  from  the  two  beam-swinging  experiments,  we 
obtain  information  on  0 ( K ) for  K vertical  as  well  as  on 
the  vertical  component  of  the  fluctuating  air  velocit \ 

tv.. 

Before  we  proceed  to  consider  the  experimental  re 
suits,  the  details  of  which  are  given  in  G jessing  and 
Bjfrresen  ( 1 ‘X> 7 ) , let  us  briefly  discuss  the  radiosonde 
experiment. 

3.  The  radiosonde  experiment 

Conventional  radiosondes  were  used  in  connection 
with  the  beam-swinging  experiments  over  a height  in- 
terval from  500  m to  some  11  km  to  obtain  the  mean 
vertical  profiles  of  temperature,  humidity,  wind  speed 
and  wind  direction. 

4.  Experimental  results 

As  shown  earlier,  the  two  experiments  give  <M  K)  for 
K vertical  and  the  vertical  component  of  the  fluctuating 
air  velocity  5K,.  Fig.  1 shows  the  variability  of  the 
spectrum  slope  n which  is  the  exponent  in  the  law 
0(A’) « A." Note  that  the  median  value  of  n for  the 
total  ensemble  is  very  close  to  — 1 1 3 in  accordance  with 
what  would  be  expected  from  theoretical  arguments 


when  the  assumption  is  made  that  we  are  dealing  with 
the  inertial  subrange  ixirtion  of  the  spectrum.  In  our 
ex|)eriment  , Kj  was  varied  from  2ir  1.5  m to  2*  0.5  m 
using  the  most  distant  receiver  (Spodsbjerg)  and  from 
— ir  3.4  m to  2*  0.7  m using  the  Utrup  receiver.  Fig.  1 
also  shows  that  in  the  height  range  from  2200-11,000 
m,  the  slope  >i  is  quite  variable  whereas  this  variability 
is  far  less  pronounced  in  the  height  interval  from  (>00- 
4700  m. 

From  meteorological  considerations  one  would  expect 
the  degree  to  which  the  atmosphere  is  mixed  to  be 
larger  in  the  lower  height  interval  where  ground  friction 
should  be  of  importance.  It  is  reasonable,  therefore,  to 
expect  the  results  from  the  spectral  measurements  to  be 
in  closer  accord  with  the  classical  11  3 law  in  this 
height  interval  than  would  the  results  from  higher  ele- 
vations where  a significant  vertical  total  stability 
generally  prevails. 

To  demonstrate  this  property  , a correlation  analysis 
was  performed  Specifically,  the  slope  »i  was  correlated 
with  various  atmospheric  parameters  such  as  the 
Richardson  number  Ri.  the  Yaisala  Brunt  frequency 
which  is  the  numerator  of  Ri.  and  the  Showalter  index 
which  is  the  differ*  nee  in  |M>tcntial  temperature  between 
the  5(H)  and  the  850  mb  surfaces.  These  atmospheric 
parameters  were  obtained  from  the  conventional  radio- 
sonde data,  the  radiosonde  being  released  so  as  to  pass 
through  the  scattering  volume  of  the  200  mi  Ktisti- 
ansand  Spodsbjerg  path  generally  at,  or  close  to,  the 
time  when  the  beam  swinging  experiments  were  per- 
formed. Since  the  radiosonde  (conventional  V’.  S. 
Weather  Bureau  type)  measures  the  pertinent  param- 
eters at  discrete  height  levels,  all  data  obtained  at  each 
height  were  fed  into  a digital  computer  XCTH'  3600) 
for  further  processing.  It  was  found  that  the  Richard- 
son number,  in  general,  was  quite  variable,  drastic 
changes  being  generally  observed  from  one  height  in- 
terval to  the  next  as  a result  of  very  small  but  variable 
vertical  wind  shears  dv  Jz  which  appear  in  the  denomi- 
nator of  the  expression  for  Ri,  i.e., 

-(—+-) 

T \ ds  (J 

Ri  — . 


Here  g is  the  gravitational  constant,  T the  mean  tem- 
perature (for  the  pertinent  height  interval)  and  cp  the 
specific  heat  at  constant  pressure. 

The  beam  swinging  experiment  measures  the  average 
value  of  the  <>(A  ) function  in  the  height  interval  from 
22(H)  11,000  m (in  the  case  of  the  Spodsbjerg  results). 
I'lius.  we  get  one  average  value  of  the  exponent  >i  for 
this  height  interval.  As  a result,  in  order  to  correlate  n 
with  appropriate  atmospheric  stability  parameters,  the 
average  value  of  each  such  parameter  must  be  com 
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I n:  1.  Histogram  of  the  distribution  of  n,  the  exponent  in  the  relation  0 1 A ' ■*  A ",  as 
obtained  at  two  receiving  stations  some  100  and  -00  mi,  respectively,  from  the  transmitter. 
I'he  heavy  line  is  the  consolidated  results  from  both  stations. 


puted  for  each  balloon  sounding.  Since  Ri— as 
Jv/dt  — * 0 and  since  a large  number  of  radiosonde 
soundings  showed  height  intervals  where  Jr  its  was 
virtually  zero,  Richardson’s  number,  as  applied,  cannot 
be  expected  to  be  a convenient  parameter  for  the  corre- 
lation analysis. 

By  using  a modified  version  of  Ri,  however,  where 
the  influence  of  (Jt  J*)*  on  Ri  was  very  much  reduced, 
the  dynamic  term  Jt  Js  did  contribute  to  the  correlation 
between  n and  the  modified  stability  parameter,  but  .is 
it  was  almost  negligible  it  was  omitted  for  the  sake  of 
simplicity  from  the  final  multiple  regression  analysis. 


I'he  analysis  thus  involved  the  new  parameter 
. £/‘ir  S \ 

,,*=—(• — H H).,)57«o  )> 

1 \ its  i p / 

" here  v,*  is  a modified  version  of  the  \ aisalit  Brunt  fre- 
quency, 7 the  average  temperature  in  the  pertinent 
height  interval  (from  2.7  11  km  for  the  Spodsbjerg  re 
suits),  JT  Jz  the  corresponding  average  temperature 
gradient,  and  7Yn  the  temperature  at  the  800  mb 
surface. 

big.  2 shows  the  result  of  the  correlation  analysis 
between  u and  < r Note  that  the  correlation  coetli 
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Kio.  2.  Result  of  a correlation  analysis  between  the  spectrum  slope  — « as  deduced  from  the 
beam-swinging  experiments  and  the  static  stability  of  the  atmosphere. 


cient  is  0.83%  for  37  ensemble  pairs.  Note  also  that  all 
the  data  obtained  at  the  Spodsbjerg  station  are  in 
eluded  except  for  cases  where  severe  radio  ducting 
saturated  the  radio  receiver,  making  a beam-swinging 
experiment  meaningless. 

Checking  the  significance  of  the  0.83  correlation 
coefficient  using  the  "student's  I test,”  one  finds  that 
the  probability  of  obtaining  a correlation  coefficient  of 
0.83  from  a sample  of  37  by  pure  chance  is  1 in  10* 
(<=  9.79).  Our  0.83  correlation  coefficient  is  thus  con- 
sidered as  highly  significant. 

The  correlation  analysis  was  repeated  on  the  basis 
of  the  Utrup  results  (height  interval  from  0.6-4. 7 km). 
For  this  ensemble  of  experimental  results,  however,  no 
significant  correlation  was  observed.  There  are  two 
obvious  practical  reasons  for  this.  From  an  inspection 
of  the  histograms  of  Fig.  1,  we  first  see  that  the  varia- 
bility (width  of  histogram)  of  the  spectrum  slope  h is 
much  smaller  for  the  Utrup  results  as  compared  with 
those  from  Spodsbjerg.  This  has  a direct  influence  on 
the  correlation  coefficient;  indeed,  a very  small  varia- 
bility of  one  of  the  parameters  involved  in  a two-pa- 
rameter correlation  analysis  makes  such  an  analysis 
meaningless.  Secondly,  in  measuring  the  vertical  pro- 
file of  temperature  through  the  scattering  volume  it  was 
observed  that  in  the  lower  part  of  the  height  interval 
involved  in  the  Utrup  experiments  the  temperature 
gradient  < IT  ds  varied  substantially  with  height.  Since 
we  require  a mean  value  of  this  gradient,  this  strong 


height  dependence  of  dT  dz  made  it  difi.cult  to  obtain 
a representative  average  profile. 

In  this  brief  presentation  of  the  results,  the  dynamical 
aspects  related  to  the  experimental  findings  will  not  be 
discussed.  We  now  consider  another  correlation  analysis 
based  on  the  beam-swinging  data  and  on  the  modified 
Yiiisala-Brunt  frequency  stability  parameter  «>i*. 

From  Eq.  (2)  we  see  that  the  slope  of  the  .Y  (speed  of 
fading  6F)  vs  6 curve  resulting  from  measurement  of  the 
Doppler  spectrum  width  SF  as  the  scattering  angle  Sis 
varied  in  such  a way  (Experiment  A)  as  to  keep  the 
direction  of  K vertical  (i.e.,>  =0),  gives  us  information 
about  the  quantity  <51,,  the  vertical  component  of  the 
fluctuating  air  velocity.  However,  since  <5U,  can  be  deter- 
mined uniquely  from  this  experiment  alone  only  when 
<5l’,  is  independent  of  height,  we  select  the  cases  where 
Experiment  B indicates  a homogeneous  atmosphere 
from  the  point  of  view  of  wind  motion,  i.e.,  where  the 
speed  of  fading  SF  is  constant  with  height.  (See  Figs. 
5 and  6 of  (i jessing  and  Byuresen.) 

The  ensemble  of  beam-swinging  results  thus  ob- 
obtained  was  used  as  the  basis  for  the  correlation 
analysis  between  the  A.Y  Afi,  which  is  proportional  to 
6l’„  and  the  stability  parameter  »•!*,  the  results  being 
given  in  Fig.  3.  We  see  that  the  correlation  coefficient 
is  0.91  and  that  large  values  of  the  vertical  fluctuating  air 
velocity  61',  are  associated  with  small  values  of  sta- 
bility, and  conversely,  that  small  values  of  61',  are 
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Fic.  3.  Result  of  a correlation  analysis  between  the  spread  in 
vertical  air  velocity  as  deduced  from  two  beam-swinging  experi- 
ments and  the  static  stability  of  the  atmosphere. 

associated  with  large  values  of  tu*.  This,  indeed,  is  to 
be  expected. 

5.  Discussion 

Having  presented  the  experimental  results  and  having 
discussed  some  of  the  important  experimental  (engi- 
neering) limitation  and  sources  of  error,  we  now  con- 
sider the  manner  in  which  the  experimental  results  are 
related  to  the  dynamics  of  the  atmosphere. 

Referring  to  Fig.  1,  we  observe  a large  variability  of 
the  refractive  index  spectrum  slope  in  the  height  in- 
terval from  2200-11,000  m.  Bearing  in  mind  that  the 
4>(K)  spectrum  is  probed  in  the  wavenumber  region 
from  2w/1.5  m to  2ir/0.5  ni  where  the  inertial  sub- 
range K~ul 5 law  is  expected  to  prevail,  this  variability 
is  somewhat  unexpected.  However,  if  we  compare 
the  current  results  with  those  obtained  using  an  entirely 
different  technique  (Eklund  el  al.  1968;  Bolgiano,  1960), 
we  find  striking  agreement  insofar  as  the  variability  of 
the  spectrum  slope  is  concerned.  Nevertheless,  in  spite 
of  this  experimental  support,  a critical  evaluation  of 
our  experimental  findings  appears  to  be  in  order.  In 
particular,  it  is  of  interest  to  know  whether  the  static 
stability  of  the  atmosphere  really  does  have  an  effect 
on  the  refractive  index  spectrum  (turbulence  spectrum) 
in  what  is  considered  as  being  the  high  wavenumber  end 
of  the  spectrum,  or  if  the  observed  correlation  (Fig.  2) 
is  merely  a result  of  the  manner  in  which  the  shape  of  the 
refractive  index  spectrum  is  deduced. 

We  have  already  observed  that  the  influence  of  the 
change  in  position  of  the  scattering  volume  has  been 
taken  into  consideration.  \..,at  remains  to  be  discussed 
is  the  theoretical  expression  with  which  the  experimen- 
tal curves  are  compared  in  order  to  obtain  a measure  of 


the  spectrum  slope  n.  If  the  refractive  index  spectrum 
can  be  written  in  the  form  0(ff)«A'_",  the  received 
power  for  a given  scattering  angle  relative  to  the  maxi- 
mum value  obtained  when  both  the  transmitting  and 
the  receiving  antennas  are  directed  horizontally  is  given 
by 

P / 2a'- 

- = (l  + — ) , 

Pt)  V 9(J  ' 

where  a is  the  elevation  angle  of  the  antennas  referred 
to  the  earth  tangent  planes  and  8t>  the  ar.'le  between 
the  tangent  planes.  If  there  are  no  vertical  mean  re- 
fractive index  gradients  causing  refraction,  then  80 =<f/ 
R,  where  R is  the  radius  of  the  earth  and  d the  path 
length. 

When  determining  the  exponent  n from  the  experi- 
mental received  power  P vs  elevation  angle  a curves, 
the  value  9„=d/R  was  used  for  all  the  experiments. 
However,  90  is  expected  to  vary  since  the  refractive 
index  gradient  dN/dz  is  expected  to  vary;  specifically, 
if  the  gradient  dN/dz  decreases,  90  increases.  Conse- 
quently, by  using  the  value  d/R  instead  of  9 o when  de 
ducing  the  exponent  | — n| , we  are  deducing  a value  of 
| — n | which  is  too  large.  The  question  now  is  whether 
this  can  account  for  the  encountered  variability  in  — « 
and  also  for  the  correlation  established  in  Fig.  2. 

Completing  our  line  of  argument,  we  see  that  if 
\dN/dz\  decreases  \dT/dz\  must  also  decrease  for  a 
given  constant  humidity,  although  the  humidity  term 
is  not  included  in  the  regression  analysis  shown  in  Fig.  2. 
Finally,  bearing  in  mind  that  the  Vaisala- Brunt  fre- 
quency v 2 is  proportional  to  the  difference  in  gradient 
between  the  adiabatic  profile  and  the  measured  tem- 
perature profile,  we  see  that  if  \dT/ds\  decreases  in  the 
stable  range,  v*  must  increase.  Consequently,  on  the 
basis  of  the  above  argument,  refraction  effects  are  ex- 
pected to  give  the  same  general  relationship  between 
e2  and  n as  that  observed  in  the  correlation  analysis  of 
Fig.  2.  As  a result  of  these  conclusions,  additional 
analyses  were  made.  The  first  involves  the  correlation 
between  the  gradient  dN/dz  as  obtained  from  the  radio- 
sonde and  the  exponent  n as  deduced  from  beam- 
swinging experiments;  no  significant  correlation  coeffi- 
cient resulted. 

Secondly,  the  minimum  scattering  angle  So  was  calcu- 
lated from  the  known  value  of  dN/dz  for  each  set  of  beam- 
swinging experiments.  Using  this  value  of  So  and  the 
experimental  power  vs  elevation  angle  relationship,  the 
exponent  « was  obtained.  Finally,  n was  correlated 
against  v?  as  in  Fig.  2.  The  results  of  this  analysis  can 
be  summarized  as  follows : 

1)  The  correlation  coefficient  was  essentially  the  same 
as  that  obtained  previously. 

2)  Points  in  Fig.  2 corresponding  to  high  values  of 
| — n\  were  shifted  slightly  down  the  — n axis  whereas 
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points  corresponding  to  low  values  of  | — n | remained 
essentially  fixed. 

In  conclusion,  therefore,  it  is  highly  probable  that  the 
slope  *i  of  refractive  index  spectrum  is  affected  by  the 
static  stability  parameter  n*. 

If  we  now  inspect  the  particular  form  of  the  current 
stability  parameter  •»>*,  we  see  that  the  temperature  at 
the  800-mb  surface  is  involved.  Excluding  this  from 
the  correlation  analysis,  the  correlation  coefficient  drops 
from  0.8.1  to  0.68;  we  thus  see  that  this  parameter  is  of 
considerable  importance.  The  quest  ion  which  now  arises 
is  whether  this  7'«oo  term  constitutes  a weighting  of 
the  density  term  f/T  such  that  the  stability  parameter 
should  be  of  the  form 

< /‘lT  K \ 


7 +c7'»(ki\ dt  Cp/ 


where  r is  a weighting  constant,  or  if  is,  in  fact,  a 
weighting  term  on  the  restoring  force  (</T/<h-h g/cp). 
This  problem  is  in  the  process  of  being  investigated. 

We  now  consider  the  correlation  analysis,  the  results 
of  which  are  presented  in  Fig.  J.  We  see  that  when  the 
contribution  to  fading  of  the  vertical  fluctuating  com- 
ponent of  the  air  velocity  is  plotted  against  si’,  a cor- 
relation coefficient  of  0.V1  results.  In  this  context  it 
should  be  born  in  mind  that  4F,,  as  deduced  from  the 
beam  swinging  experiments,  measures  the  variability 
in  the  vertical  speed  of  the  scattering  elements  at  essen 
tially  one  given  instant.  What  Fig.  .1  tells  us,  therefore, 
is  that  the  velocity  distribution  of  the  scattering  ele- 
ments becomes  narrower  as  the  static  stability  param- 
eter increases  (Gjcssing,  1062). 

This  indeed  is  to  be  expected ; if  the  stability  is 
sufficiently  high  one  would  expect  the  scattering  ele 
ments  in  a limited  spatial  region  to  oscillate  vertically 
in  unison,  the  frequency  of  oscillation  being  the  Viiisiila- 
Brunt  frequency  Under  such  conditions  one  would 
expect  4V,  to  approach  zero  as  the  spatial  region  re- 
sponsible for  the  scattering  is  limited. 

The  fact  'hat  we  are  indecal  obtaining  a clear  correla- 
tion between  the  modified  YiiisiUii-Brunt  frequency  i'i* 
and  4V,  adds  substantially  to  the  confidence  level  of 
the  — « vs  *|*  correlation  analysis  of  Fig.  2.  Our  confi- 
dence in  fi1  as  a stability  parameter,  and  not  merely  as  a 
refraction  parameter,  is  improved. 

Accepting  then,  on  the  basis  of  the  above  arguments, 
that  we  have  established  a close  relationship  between 
the  spectrum  slope  and  the  atmospheric  stability,  the 
question  which  arises  is  how  it  is  possible  for  the  static 
stability  to  affect  the  spectrum  slo|>e  in  a wavenumber 
region  which  is  expected  to  fall  within  the  inertial 
subrange. 

Several  authors  (Bean  ft  al.,  Kaimal  and  Haugen, 
1%7;  Panofsky  and  McCormick,  l%();  Panofsky, 
1^68)  have  discussed  the  shape  of  the  low  wavenumber 
end  of  the  spectrum  in  relation  to  atmospheric  stability 


and  in  relation  to  height  above  the  ground.  'Ilirse  con- 
tributions, however,  do  not  appear  to  have  any  bearing 
on  the  high  wavenumber  end  of  the  spectrum  such  as 
we  arc  currently  investigating  at  large  distances  above 
the  ground.  Based  on  a multi-frequency  radio  experi- 
ment, however,  Bolgiano  (IdoOl  has  observed  a varia- 
bility of  the  n exponent  (from  t>/i  to  18/.I)  very 
similar  to  that  encountered  in  the  current  experiment. 
Furthermore,  in  correlating  the  exponent  ti  with  a dy- 
namic stability  parameter  (Richardson's  number),  he 
obtains  a correlation  coefficient  of  0.80  from  a 20  point 
sample  and  concludes  that  there  seems  to  be  a strong 
connection  between  the  viscous  cutoff  scale  and  the 
dynamic  stability  of  the  atmosphere. 

In  this  contribution  we  shall  refrain  from  offering  an 
explanation  pertaining  to  the  relationship  established 
between  spectrum  slojre  and  atmospheric  stability;  it  is 
the  hope  that  this,  in  due  course,  will  be  given. 

0.  Conclusions 

These  conclusions  can  be  summarized  as  follows. 
Expressing  the  static  stability  of  the  atmosphere  by  a 
modified  version  of  the  familiar  Vhisidii-Brunt  frequency 
<'*,  it  is  found  that  both  the  slo|te  of  the  refractive  index 
spectrum  and  the  velocity  spread  of  the  irregularities 
(scatterers)  causing  scattering  of  radio  waves  arc 
strongly  influenced  by  the  static  stability.  The  fact  that 
the  width  of  the  velocity  distribution  (as  determined 
from  a radio  experiment)  is  so  well  correlated  with  our 
stability  parameter  (correlation  coefficient  0.11)  is 
comforting.  As  has  been  shown,  we  have  laboriously 
examined  a variety  of  phenomena  that  conceivably 
could  influence  our  interpretation  of  the  experimental 
results.  On  the  basis  of  these  examinations  the  con- 
clusion can  be  drawn  that  the  stability  is  indeed  proba- 
bly influencing  tbe  shape  of  tbe  refractive  index  spec- 
trum. and  hence  the  turbulence  field.  It  is  unlikely, 
that  the  variability  of  the  spectrum  slope  is  determined 
by  refraction  effects.  The  confidence  in  this  statement 
is  very  much  increased  by  the  fact  that  there  is  a strong 
correlation  between  stability,  as  defined,  and  the  degree 
of  irregularity  of  the  air  motion. 

What  remains  then  is  to  introduce  arguments  based 
on  turbulence  theory  with  the  aim  of  explaining  tbe 
relationships  established  in  this  paper;  this,  however, 
will  be  (he  topic  of  future  work. 
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Abstract.  Measurements  have  been  made  with  fast-response  multi-channel  temperature,  humidity 
and  refractive  index  sensors  flown  to  2000  m on  a tethered  balloon  to  investigate  small-scale  fluctu- 
ations important  in  radio-wave  scattering,  their  relation  to  atmospheric  parameters,  and  their  spatial 
variation  in  both  one  and  three  dimensions.  Data  from  the  three  types  of  sensors  at  one  point  were 
consistent  for  frequencies  up  to  about  8 Hz.  Power  spectra  of  data  at  various  heights  were  computed 
over  0.1  to  10  Hz  and  generally  showed  slopes  (on  a log-log  plot)  close  to  -5/3  above  I Hz  but  ranged 
from  — 1.5  to  —3.5  at  lower  frequencies;  in  this  range  (/<  1 Hz)  slopes  were  close  to  —5/3  for 
negative  Richardson  number  (Ri),  provided  temper;.ture  gradients  were  steeper  than  - l.l°C  100  m 1 
and  wind  shears  > 1.4  x I0  2 s 1 approx.  Steeper  slopes  were  generally  associated  with  stable  atmo- 
spheric conditions  but  no  precise  relation  to  the  above  parameters  was  found.  Spectral  density  was  a 
maximum  for  Ri  ~ —0.75. 

Cross-correlations  of  0.5  were  frequently  observed  between  sensors  1 m apart  in  orthogonal 
directions;  in  the  vertical,  examples  of  negative  correlation  of  vapour  pressure  were  occasionally 
found  over  this  spacing.  Using  four  sensors  spaced  in  line  over  9 m,  cross-spectrum  phase  calcu- 
lations of  drift  speeds  were  found  to  be  consistent  with  measured  wind  speeds.  The  ratio  of  ‘iden- 
tification distance'  (coherence  = 0.6)  to  scale  size  of  irregularities  ranged  from  0.25  to  0.5  with  no 
apparent  relation  to  height  or  meteorological  parameters. 


1.  Introduction 

This  paper  describes  measurements  of  temperature,  humidity  and  refractive  index 
fluctuations  in  the  troposphere,  using  fast-response  instruments  flown  on  a tethered 
balloon,  which  have  been  made  with  the  object  of  studying  small-scale  structure  which 
affects  radiowave  propagation.  The  scattering  of  radio  waves  in  the  troposphere 
depends  on  the  variation,  in  space  and  time,  of  the  radio  refractive  index;  the  presence 
of  small-scale  variations  of  refractive  index,  important  in  this  scattering  process, 
results  from  small-scale  fluctuations  of  temperature  and  humidity. 

Radio  scattering  experiments  designed  to  measure  the  three-dimensional  spatial 
spectrum  of  refractive  index  variations  in  the  troposphere  show  that  the  spectrum 
varies  substantially  with  the  dynamic  stability  of  the  troposphere  (Gjessing  el  al., 
1969).  It  is  found  that  whereas  the  classical  Obukhov-Kolmogorov  theory  predicts  a 
wave-number  power  spectrum  in  the  inertial  subrange  of  the  form  E(k)~k~"  with 
n~ 5/3,  radio  experiments,  and  indeed  direct  measurements  in  the  troposphere,  show 
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that  n may  vary  from  I to  about  3.  Moreover,  it  is  still  not  clear  to  what  extent  the 
fluctuation  characteristics  (scale  si/e,  spectral  slope)  of  one  parameter  for  example, 
humidity  can  be  predicted  from  measurements  of  another  parameter.  A knowledge 
of  the  instantaneous  small-scale  spatial  structure  is  also  important.  To  study  these 
features,  a spaced-probe  array  consisting  of  a multi-channel  refractometer  with 
temperature  and  humidity  sensors  attached  was  flown  on  a tethered  balloon  at  heights 
up  to  2000  m.  The  spatial  arrangement  of  the  sensors  could  be  varied,  and  different 
patterns  were  used  to  suit  particular  objectives.  These  were  primarily  to  investigate 
the  relation  between  refractive  index,  humidity  and  temperature  fluctuations  at  one 
point;  to  measure  spatial  variations  in  three  dimensions,  and  to  measure  the  cor- 
relation between  sensors  spaced  in  one  horizontal  line.  The  measurements  were  made 
during  July  and  August  l%4  in  anticyclonic  conditions,  giving  stable  and  near- 
neutral  lapse  rates  and  low  wind  speeds.  Under  these  conditions,  fluctuations  were 
often  small  and  intermittent.  In  the  following  sections,  significant  and  unusual  results 
from  an  analysis  of  the  various  data  are  described 

2.  Fxperinicntal  Procedure 

The  experiments  were  carried  out  using  the  tethered  balloon  lacilitv  at  ( aldington. 
Bedfordshire.  ILK.  where  it  was  possible  to  fly  instruments  at  heights  up  to  '*>00  in 
The  equipment  consisted  of  an  \-band  microwave  refractometer.  line  wire  therm*' 
meters,  fast-response  humidity  sensors,  a cup  anemometer  and  a radio  lelemetrv 
system.  The  refractometer  had  four  open  sampling  cavities  and  a closed  reference 
cavity  in  the  swept-frequency  Hirnbaum  system  ( Hirnhaiim.  losoi  ,,i„|  was  developed 
from  an  earlier  model  described  by  fowler  <7  ul  ( I'X'M  It  has  a Irequencv  response 
of  0 to  20  11/  (3-dB  point).  The  temperature  sensors  were  made  from  10 n platinum 
wire  and  were  used  in  a resistance-thermometer  circuit,  giving  a tesponse  to  so  11/ 
The  humidity  sensors  were  of  the  coated  quartz-crystal  tv  pc  capable  ol  giving  Iasi 
response,  to  25  Hz  approximately.  Both  the  temperature  and  hunuditv  sensors  have 
been  described  in  detail  by  Gjcssing  et  ol  (l')r>S) 

(file  each  of  the  thermometers  and  humidity  sensors  was  placed  close  to  a re 
fractometer  cavity  so  that  consistent  fluctuation  measurements  might  be  expected  up 
to  some  20  11/  at  normal  wind  speeds  Three  different  arrangements  of  the  sensor 
groups  were  used;  these  were,  firstly,  an  orthogonal  array  with  one  group  ol  sensors 
at  the  origin  and  one  on  each  axis  at  equal  distances  of  I m.  secondly,  the  sensor 
groups  were  placed  in  line  on  a *f-m  boom,  flown  so  as  to  be  in  line  with  the  mean 
wind  as  determined  by  a wind  sock  attached  to  the  end  of  the  boom  The  spacings 
were  in  order  1.5,  4.4,  3.2  m,  which  gave  six  different  spacings  between  pairs  for 
spatial  correlation  measurements.  The  third  arrangement  had  sensors  placed  on  a 
4.5-m  boom  cross-wind,  with  spacings  of  0.75,  2.3  and  1.5  m 
The  data  were  telemetered  to  the  ground  using  a standard  IRK)  multiplex  FM 
subcarrier  system  modulating  a 1500-MH/  FM  transmitter.  The  multiplex  was 
recorded  on  magnetic  tape  and  analog  signals  were  recorded  on  paper  charts 
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The  data  were  analysed  by  standard  digital  computer  techniques;  selected  data 
were  digitised  at  40  samples  per  second  after  the  signals  had  been  filtered  to  the 
required  bandwidth;  using  the  methods  of  Blackman  and  Tukey  ( I *>5S>  the  auto-  and 
cross-correlations  of  the  signals  were  calculated,  together  with  power  spectra  ,V(  /') 
and  cross-spectral  properties.  These  included  co-spectra  (P, :)  and  quadrature  spectra 
(£>,.),  coherence  t'OII,  , = t Qt:)/S ,S2,  i.e.,  the  normalised  magnitude  squared 
of  the  cross-spectrum  between  channels  I and  2,  and  the  associated  phase  shift 
(*, . tan  Prewhitening  and  recolouring  processes  were  used  in  computing 

spectra.  Digital  filtering  techniques  were  also  used  to  remove  dc  and  low  frequency 
components  which  can  give  rise  to  errors  in  analysis,  and  to  eliminate  high  frequencies 
and  reduce  the  number  of  data  points  when  studying  low-frequency  components. 

3.  Sensor  Performance 

The  first  point  of  interest  is  to  compare  the  fluctuations  measured  simultaneously 
by  each  type  of  sensor  at  one  location  and  to  see  if  they  are  consistent,  figure  I shows 
such  fluctuations  at  each  of  four  locations  and  illustrates  the  mutual  consistency  of 
the  recorded  data.  These  were  recorded  at  1500-m  height  with  the  9-nt  in-line  array 
and  show  small  fluctuations  and  occasional  large  changes,  intermittcncy  typical  of 
stable  conditions;  the  cross-correlations  between  sensors  of  the  same  type  exceeded 
0.95  for  all  channels. 

In  f igure  2a  small-scale  measured  refractive-index  fluctuations  typical  of  more 
uniform  turbulence  near  the  ground  are  compared  with  those  calculated  from  the 
humidity  and  temperature  fluctuations;  the  agreement  is  clear  and  the  cross-cor- 
relation is  0.42  over  a 2-ntin  sample.  The  correlation  is  thought  to  he  limited  mainly 
by  non-linearity  of  the  humidity  sensors  particularly  at  large  relative  humidity 
(Gjessing  <■;  a I .,  1468).  The  calculated  values  were  obtained  using  the  equation 

SN  = aSe  -f  fiST , (I) 

where  <V(mb)  is  the  fluctuation  in  water  vapour  pressure  and  d 7 ( C)  the  fluctuation 
in  temperature.  This  equation  is  obtained  bv  differentiation  of  the  Debye  equation 
in  the  reduced  form 


where  N is  one  million  times  the  amount  by  which  the  radio  refractive  index  exceeds 
unity,  /’(mb)  is  the  total  pressure,  7 (K)  is  the  temperature  and  < (mb)  is  the  water 
vapour  pressure  (e  g.,  Bean  and  Dutton.  1966).  Typically  a - 4.5  mb  and  /i  - —1.2  t . 

f igure  2b  shows  the  power  spectra  to  be  similar  to  at  least  b If/.  The  root  coherence, 
representing  the  cross-correlation  in  each  frequency  band  of  the  spectrum,  is  shown 
in  Figure  2c  and  is  greater  than  0.7  up  to  some  o Hr.  The  dip  at  7 If;  is  thought  to 
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Fig.  I.  Examples  of  refractivity  N,  relative  humidity  H and  temperature  7"  fluctuation,  9-m  in-line 
array,  1500  m,  1969  July  30  at  1448. 


be  caused  by  a weak  7-Hz  noise  component  present  in  the  refractometer  output. 
Above  this  frequency,  the  decrease  in  coherence  may  be  due  to  fluctuations  in  wind 
direction,  and  to  finite  sensor  spacings  of  about  20  cm.  However,  the  general  agree- 
ment between  sensor  measurements  is  considered  acceptable. 

4.  Single-Point  Measurements 

Before  considering  the  results  from  spaced-array  measurements,  much  information 
can  be  deduced  from  the  observations  at  one  point  in  relation  to  atmospheric  structure. 
Samples  of  atmospheric  fluctuations  were  over  periods  of  5 to  10  min  at  heights  up 
to  2000  m,  often  at  150-m  intervals,  and  from  these  we  have  sought  answers  to  our 
particular  questions: 
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Fig.  2.  (a)  Comparison  between  measured  refractivity  fluctuation  and  that  calculated  from  tem- 
perature and  humidity  fluctuations  at  one  point.  Test  data  2 m above  ground.  - (b)  Corresponding 
power  spectra.  - (c)  Coherence  between  measured  and  calculated  fluctuations. 

(a)  To  what  extent  do  the  spectra  S(f)  of  the  N,  T and  c fluctuations,  as  measured 
directly,  vary  with  time  and  position? 

(h)  To  what  extent  do  the  slopes  n (on  a log-log  plot)  of  spectra  S(f) at/-"  vary 
in  unison? 

(c)  To  what  extent  do  the  standard  deviations  <r,,  <t,,  <tv  correlate? 

(d)  To  what  extent  do  n and  a correlate  with  conventional  parameters  describing 
the  stability  of  the  atmosphere,  such  as  the  Richardson  number  Ri.  temperature 
gradient  dT/dz,  and  wind  shear  dK/dz? 

4.1.  Variation  of  spectra  in  space  ano  time 

Typical  e and  T power  spectra  (c.g..  Figure  8)  show  that  the  high-frequency  end  of 
the  spectrum  from  I to  10  Hz  nearly  always  follows  a -5/3  power  law  irrespective 
of  the  position  in  space  at  which  the  measurements  are  made.  This  also  applies  to 
the  N spectra  since  the  fluctuations  are  linearly  related.  However,  the  N,  e,  T spectrum 
slopes  are  found  to  vary  considerably  for  frequencies  below  some  1 Hz.,  usually  being 
steeper  than  —5/3,  and  the  frequency  at  which  the  spectrum  slope  changes  varies 
with  time  and  position  in  space.  When  reduced  to  scale  size  using  the  mean  wind 
speed,  this  point  varies  over  I to  10  m.  The  spectrum  slopes  arc  also  found  to  vary 
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considerably  with  time  and  position  in  space.  This  is  probably  a consequence  of 
conditions  prevailing  during  most  of  the  experiments,  namely  fine  weather  often 
associated  with  strong  convection,  low  wind  speed  und  some  degree  of  layering  as 
shown  by  radiosonde  ascents.  In  these  conditions,  intermittency  was  frequently  a 
marked  feature  of  the  data. 

4.2.  Variation  of  spectral  slopes 

In  answer  to  question  (b),  Figure  3 shows  the  slopes  of  measured  humidity  and  re- 
fractive index  spectra  plotted  against  temperature  spectrum  slopes  for  frequencies 
in  the  range  0. 1 to  I Hz.  The  slopes  vary  over  a wide  range,  but  the  plotted  points 
tend  to  lie  along  a 45'-  line  indicating  that  on  many  occasions  the  spectrum  slopes  are 
well  represented  by  that  of  the  corresponding  temperature  spectrum  alone.  There  is 
no  systematic  height  dependence  of  the  data  in  Figure  3. 


Fig.  3.  Slopes  of  measured  humidity  and  refractive  index  spectra  against  temperature  spectrum 
slope  in  the  range  0.1  to  1.0  Hz.  Heights:  0 to  1500  m. 


4.3.  Correlation  of  standard  deviations 

There  was  some  evidence  of  a weak  correlation  between  the  standard  deviations 
<rr,  a,  but  instrumental  shortcomings  prevented  more  positive  conclusions  from 
being  reached.  This  arises  as  follows.  The  contribution  to  a is  largely  from  the  low 
frequency  end  of  the  spectrum  where  the  spectral  density  is  large.  Therefore  in  order 
to  get  an  accurate  estimate  of  a„  one  must  know  the  calibration  curve  for  the  humidity 
sensor  very  accurately  over  the  whole  range  of  humidities.  Since  this  curve  is  non- 
linear, with  the  slope  increasing  substantially  with  relative  humidity  (Giessing  et  a /., 
1968),  small  variations  in  the  calibration  during  an  experiment  will  have  a large  effect 
on  the  o,-o T correlation. 

It  is  also  of  interest  to  examine  the  standard  deviations  of  the  three  measured 
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parameters  under  conditions  when  the  temperature  and  humidity  contributions  to 
refractive  index  fluctuations  are  of  similar  magnitude.  On  many  occasions,  partic- 
ularly on  the  large  scale,  the  humidity  fluctuations  dominate  but  some  examples  were 
obtained  where  the  small-scale  fluctuations,  equivalent  to  frequencies  greater  than 

0 I H/,  fulfilled  the  former  requirement.  While  these  are  related  linearly  (Equation  (I )) 
their  instantaneous  relative  magnitudes  depend  upon  physical  processes  causing  the 
fluctuations.  If  the  parameters  are  considered  as  independent  and  passive,  carried  in 
the  turbulent  field,  then  a,  and  a,  contribute  randomly  to  <rN;  but  if  the  fluctuations 
are  related,  the  variances  are  correlated,  and  from  Equation  (I)  we  derive 

CTv  = + 23pCrT , (31 

where  C,r  is  the  covariance;  similarly  in  deriving  spectra,  the  cospectrum  must  be 
included  as  has  been  done  by  Ciossard  (I960).  Of  the  samples  studied,  generally  ot 
2-  to  5-mtn  duration,  in  no  case  did  T and  e fluctuate  in  such  a way  as  to  cancel  ot 
greatly  reduce  their  combined  contribution  to  the  variance  of  ,V.  However  this  pos- 
sibility in  atmospheric  data  is  not  ruled  out  and  can  certainly  occur  in  surface  mea- 
surements, as  indicated  by  the  data  of  Swinbank  and  Dyer  (1967)  obtained  under 
conditions  of  strong  surface  evaporation. 

4.4.  Relation  of  spectral  slope  ano  variance  to  stability  parameters 

To  investigate  the  relation  between  V,  e,  7'  and  the  dynamic  state  of  the  troposphere, 
the  Richardson  number  Ri,  wind  shear  dl  ,d;  and  temperature  gradient  d T/d:  were 
first  evaluated  using  data  from  routine  Balthum  meteorological  soundings  and 
radiosondes  which  gave  mean  temperature  and  wind  at  150-m  intervals  near  the  time 
of  the  observations.  Then  the  average  was  taken  of  the  parameters  in  the  two  height 
intervals  above  and  below  the  point  of  observation  before  an  attempt  was  made  to 
relate  these  rather  large-scale  quantities  to  the  spectral  measurements.  Of  these  data 
the  wind  shear  was  the  most  uncertain,  but  encouragement  was  given  by  some 
similarity  between  wind  profiles  at  6-h  intervals,  and  the  emergence  of  significant 
correlations  parameters  as  shown  below. 

The  relation  between  wave-number  spectrum  5(A  ) for  A = 2>r/3m~ 1 and  Ri  is  shown 

in  Figure  4.  The  spectral  density  peaks  at  Ri 0.75  both  for  e and  T and  the  curve 

is  similar  in  shape  to  that  found  by  Bean  cl  al.  (1967)  obtained  at  larger  wave  numbers 
using  tower  data  and  attributed  by  them  to  convective  mixing.  For  large  positive  Ri,  the 
turbulence  is  damped  by  static  stability  and  as  Ri  =0  is  approached,  mixing  in  a region 
of  small  potential  gradient  produces  little  enhancement  of  scalar  fluctuations.  For 
increasing  negative  Ri,  fluctuation  intensity  increases  and  the  effect  of  thermals  may 
be  important.  Note  that  in  Figure  4 the  peak  corresponds  to  measurements  taken  at 
8-m  height  and  that  the  mean  lapse  rate  from  ground  to  15  m was  2.8  V'C  100  m ' 1 ; 
however,  there  is  no  systematic  height  variation  of  the  plotted  points. 

Plotting  the  slope  n of  the  icmperaturc  spectrum  at  the  lower  frequencies,  0.1  to 

1 Hz,  against  Ri  as  illustrated  in  Figure  5a.  shows  that  n remains  close  to  5/3  for 
negative  Ri.  When  Ri  becomes  positive,  n increases  substantially,  varying  somewhat 
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Fig.  4.  Variation  of  spectral  density  of  temperature  and  humidity  fluctuations  with  Richardson 
number  Ri  at  k = 2n  3 m-1.  Heights  are  shown  in  metres  for  larger  densities.  Smaller  densities  were 
measured  between  30-  and  I500-m  height. 

unsystematically  between  2 and  3.5.  It  can  be  inferred  therefore  that  large  negative 
exponents  in  the  spectrum  power  law  are  associated  with  stable  atmospheric  con- 
ditions, but  no  precise  relationship  between  n and  Ri  can  be  deduced  for  such  cases. 
Steep  slopes  were  also  found  by  Lane  and  Paltridge  (1968,  Figure  6)  under  low 
wind-speed  conditions. 

If  instead,  n is  plotted  against  dT/dz,  again  for  / <1  Hz  as  shown  in  Figure  5b,  it 
is  found  that  n lies  close  to  5/3  for  temperature  gradients  that  are  'deeper  than  - 1 . 1 °C 
100  m *,  that  is,  slightly  steeper  than  the  adiabatic  lapse  rate.  For  smaller  gradients,  n 
varies  unsystematically  again  between  2 and  3.5.  The  relation  of  n to  wind  shear 
follows  a similar  pattern,  withn~5/3  for  shears  greater  than  1.4  x 10~2  s'1  and  with  n 
variable  under  conditions  of  weak  shear,  in  broad  agreement  with  the  conclusions 
of  Bardeau  and  Mism6  (1969)  using  radiosonde  observations.  These  variations  are 
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Fig.  5.  Variation  of  slope  of  temperature  spectrum  (/<  1 Hz)  with  (a)  Richardson  number,  (b) 
temperature  gradient  and  (c)  wind  shear. 

generally  consistent  with  the  behaviour  of  n with  Ri  and  we  may  conclude  that,  to  the 
extent  that  values  of  Ri,  dT/dz,  dV/dz  averaged  over  150  m are  descriptive  of  the 
dynamic  stability  in  the  region  of  observation,  they  give  an  indication  of  the  con- 
ditions under  which  -5/3  law  turbulence  exists. 

An  investigation  of  the  behaviour  of  standard  deviation  and  spectrum  slope  has 
failed  to  show  any  clear  relation  between  the  two,  and  this  implies  that  spectral 
density  plots  S(f)  do  not  in  general  converge  at  either  the  low  or  high  frequency  ends 
of  the  spectrum. 


5.  Results  from  the  Three-Dimensional  Orthogonal  Array 

Time-lagged  cross-correlations  were  calculated  between  pairs  of  sensors  in  the  1-m 
spacing  array.  In  general  the  maximum  correlations  covered  a wide  range  at  different 
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heights  and  times  in  the  atmosphere,  but  rarely  exceeded  0.8;  frequently  they  were 
less  than  0.5  and  were  sometimes  negative,  indicating  in  general  that  the  scale  size  of 
irregularity  structures  making  a major  contribution  to  the  variance  was  often  of  the 
order  of  a few  metres.  Compared  with  the  maximum  cross-correlation  in  the  vertical 
and  cross-wind  directions,  that  for  the  pair  of  sensors  in  line  with  the  mean  wind  was 
generally  larger;  this  is  to  be  expected  and  indicates  that  the  irregularity  structure 
does  not  change  very  rapidly  within  a travel  distance  of  I m. 

Figure  6 shows  a more  unusual  occurrence  which  is  discussed  because  of  its  interest 
and  importance  to  radio  meteorological  considerations.  Data  from  the  vertical  pair 
of  humidity  sensors  are  negatively  correlated,  possibly  representing  a situation  where 
one  of  the  sensors  penetrated  a sharp  gradient  in  vapour  pressure.  Inspection  of  the 
original  charts  showed  that  the  relative  humidity  variations  were  large-scale  and 
generally  in  antiphase  whereas  the  temperature  variations  were  faster  and  generally  in 
phase,  as  were  the  refractive  index  fluctuations.  These  observations  could  have  been 
caused  by  a thin  undulating  layer,  such  that  on  the  average  when  one  of  the  two 


Fig.  6.  Lagged  cross  correlations  between  pairs  of  sensors  in  a three-dimensional  array,  1 -m  spacing, 

900  m,  1969  July  18  at  1443. 
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vertically  spaced  sensors  was  inside  the  layer,  the  other  was  outside;  it  would  be  very 
desirable  to  have  more  information  on  the  structure  of  the  atmosphere  to  help  explain 
these  observations.  However,  of  the  available  data,  the  mean  wind  shear  and  the 
temperature  gradient  averaged  over  150-m  intervals  (used  earlier)  are  of  too  large  a 
scale  to  be  very  useful,  while  the  records  obtained  from  the  balloon  ascent  with  the 
equipment  contain  mixed  height  and  time  variations,  apart  from  changes  in  the 
height  of  the  fine  atmospheric  layer  structure  with  time.  The  only  general  conclusion 
reached  is  that  all  cases  of  negative  vertical  vapour-pressure  correlation  occurred 
when  the  larger  scale  (150  m)  parameter  indicated  stable  atmospheric  conditions.  To 
obtain  adequate  information  in  the  future  to  help  explain  this  type  of  phenomenon, 
it  would  seem  necessary  to  make  measurements  of  wind  and  temperature  at  least  at  a 
number  of  points  spaced  vertically  about  the  main  instrumentation,  perhaps  spaced 
over  some  50  m. 

Figure  7 shows  maximum  cross-correlations  between  pairs  of  sensors  in  the  three- 
dimensional  array  at  various  heights  during  one  ascent.  Two  height  regions  where  the 
vertical  humidity  sensors  show  negative  correlation  are  evident,  and  although  the 
various  correlations  cover  a wide  range,  there  is  some  systematic  height  variation  in 
the  data.  If  we  apply  Equation  (1)  to  the  sensors  at  each  of  two  spaced  points  1 and  2, 
and  assume  the  coefficients  a,  p apply  at  both  points,  then  the  covariance  C*  between 
refractometer  cavities  is  given  by: 

C*  - ot 2C.  + P1CT  + aP(C,lT1  + C,2T1),  (4) 

where  the  CeT  are  now  crossed  covariances  between  the  humidity  and  temperature 
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Maximum  cross-correlations  between  pain  of  sensors  in  a three-dimensional  array, 
1-m  spacing,  at  various  heights.  1969  July  18. 


Fig.  7. 
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sensors.  Replacing  covariances  by  cross-correlations  e.g.,  C,  = a, : CC,,  etc.,  we  obtain 

^•CCv  = Cv  = a V,CC,  + p:o2rCCT  + zpo,oT(CC,XTi  + CC,jri)  (5) 

This  equation  indicates  that,  within  bounds  determined  by  the  relevant  standard 
deviations,  the  correlations  between  sensor  pairs  can  vary  widely.  In  particular,  for 
the  data  of  Figure  6 the  positive  CCS  and  negative  CC,  are  consistent  for  a range  of 
negative  CC,T  cross-correlations,  but  as  the  coefficients  in  Equation  (5)  involve  the 
sensor  calibrations  which  it  has  been  explained  are  somewhat  uncertain,  a more 
precise  comparison  is  not  possible. 

A study  of  the  time  lags  for  maximum  cross-correlations  between  sensors  indicates 
that  the  direction  of  general  air  movement  through  the  array  varied  widely,  often  over 
± 45“  both  in  the  vertical  and  horizontal,  expccially  under  low  wind  conditions 
(1-5  m s~l).  In  particular,  vertical  velocity  components  of  4 to  I ms-1  over  periods 
of  minutes  both  upward  and  downward  on  different  occasions  were  observed.  It  is 
thought  that  these  movements  are  associated  with  rising  thermal  structure  and  the 
order  of  measured  velocity  agrees  well  with  measurements  of  the  ascent  rate  of  large 
convective  cells  by  sensitive  radar  methods  (Ottersten,  1969).  These  movements  are  of 
some  consequence  in  that  it  is  common  practice  to  make  use  of  a mean  wind  speed  r 
to  convert  frequency  spectra  S(  f)  into  wave-number  spectra  S(k)  by  k = 2n  t /T;  a 
steady  deviation  of  wind  direction  from  that  assumed  could  lead  to  errors  of  up  to 
50°o  in  k.  Also  because  of  these  movements  a simple  interpretation  of  the  correlations 
in  terms  of  characteristic  blob  shape  and  size  does  not  seem  easy  and  other  inter- 
pretations are  being  sought. 

6.  Results  from  the  Cross-Wind  One-Dimensional  Array 

In  this  experiment  a 44-m  boom  was  flown  normal  to  the  mean  wind  by  means  of 
a specially  arranged  small  parachute,  and  with  four  sensors  of  each  type  unequally 
spaced,  as  described  in  Section  2,  a total  of  6 different  spacings  between  pairs  was 
possible.  From  the  data,  lagged  cross-correlations  were  calculated  and  Figure  8 shows 
examples  of  the  spatial  correlograms  obtained  using  the  maximum  cross-correlations 
for  each  spacing;  these  maxima  do  not  always  occur  at  zero  lag,  due  to  fluctuations 
in  wind  direction.  On  the  left  (8  m),  the  correlogram  falls  rapidly,  indicating  a scale 
size  (CC=J)  of  order  14  m near  the  ground,  similar  to  the  ground-level  results  of 
Thompson  er  al.  (1968);  the  example  to  the  right  (150  m)  is  more  typical  of  obser- 
vations at  greater  heights  and  shows  a predominance  of  larger  scale  sizes. 

Also  shown  on  each  plot  is  the  average  of  the  4 autocorrelation  curves  obtained 
from  the  four  sensors  and  it  is  of  interest  to  compare  the  spatial  separations  L with 
the  time  lags  t which  give  equal  correlation,  to  obtain  a speed  V = Lit.  If  the  struc- 
ture drifting  past  the  sensors  is  isotropic  and  not  changing  with  time,  V will  be  the 
true  drift  speed,  but  if  the  structure  is  changing,  the  calculated  speed  will  be  an  over- 
estimate of  the  true  velocity.  Doing  this  computation  for  each  point  on  the  spatial 
correlograms,  the  speeds  shown  at  the  top  of  Figure  8 were  obtained. 


t 


• io  f l i to  r 

rflfQUCNCV  W HttTI  fRIOUtNCT  IN  HIRTI 


Fig.  8.  (a)  Typical  examples  of  humidity  and  temperature  spectra.  - (b)  Temperature  auto  and 
spatial  correlograms.  Cross-wind  array  1969  August  I.  (c)  Apparent  wind  speed  derived  from  (b) 
using  V Lit  where  R(l.)  RU).  - The  left-hand  curves  correspond  to  8-m  height  and  the  right- 

hand  curves  to  150-m  height. 

At  8 m above  the  ground  for  correlations  greater  than  0.2,  the  speeds  obtained  are 
nearly  all  equal  with  a mean  of  0.75  ms1  while  records  of  the  wind  speed,  measured 
with  a cup  anemometer  showed  a large  variation  ranging  from  I to  5ms'1  with  a 
mean  of  l-J  m s'1.  Thus  the  agreement  is  poor,  tiie  correlation  method  under-estimat- 
ing the  mean  speed,  and  could  be  caused  by  anisotropy  or  by  a vertical  wind.  At  150  m 
calculated  speeds  are  more  scattered  but  with  several  about  2ms'1;  the  measured 
speed  again  varied  greatly,  from  I J to  5 m s'1  with  a mean  of  2$  m s ' 1 giving  more 
satisfactory  agreement.  It  is  clear  that  conditions  of  small  and  very  variable  wind 
speed  are  unsatisfactory  for  attempting  measurements  of  spatial  structure  using 
spaced  arrays  of  sensors,  but  unfortunately  these  conditions  occurred  frequently  at  the 
time  of  our  field  experiments. 

Also  shown  in  the  lower  part  of  Figure  8 are  the  corresponding  spectra  which  were 
mentioned  earlier.  Of  some  24  spectra  from  this  series  of  measurements,  all  but  3 
exhibited  a — 5/3  power  law  at  frequencies  above  approximately  I Hz,  and  at  lower 
frequencies  the  slope  was  variable,  except  when  measured  close  to  the  ground  where  it 
was  again  -5/3  (e  g..  Figure  8);  the  steeper  slopes  were  in  all  cases  associated  with 
large  spatial  correlation  distances.  This  conclusion  was  also  reached  in  the  case  of  the 
three  dimensional  array  measurements.  There  was  some  evidence  that  the  correlation 
Its'  "'xc>  of  water  vapour  pressure  were  smaller  than  those  of  temperature. 
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7.  Results  front  the  In-Line  One-Dimensional  Array 

in  this  experiment  a 9-m  array  of  sensors  was  flown  with  the  object  of  studying  the 
changes  and  apparent  speed  of  irregularities  drifting  past  them  in  the  mean  wind. 
Again,  four  sensor  groups  were  used  with  six  different  spacings  between  pairs  as 
described  in  Section  2.  Difficulty  was  experienced  with  the  array  boom  oscillating 
under  low  wind-speed  conditions  but  some  useful  data  were  acquired.  Figure  9 shows  a 
typical  spatial  correlogram  R(L)  plotted  from  the  maximum  cross-correlation  at  each 
spacing.  Also  shown  is  the  root  coherence,  or  cross-correlation  in  each  frequency  band, 
obtained  by  applying  a Fourier  transform  to  the  lagged  cross-correlation  function  of 
each  pair  of  sensors;  this  indicates  how  the  various  frequency  components  in  the 
irregularity  spectrum  contribute  to  It  is  the  low-frequency  components  which 
make  a major  contribution:  at  higher  frequencies  the  coherence  decreases  rapidly 
with  increasing  spacing.  Using  these  curves,  it  is  of  interest  to  examine  the  relation 
between  frequency  (/)  and  the  distance  at  which  the  root  coherence  falls  to  say  0.6. 
This  has  been  done  for  data  at  different  heights  in  the  atmosphere  and  Figure  10 
shows  that  an  approximately  linear  relation  exists  between  this  ‘identification'  dis- 
tance and  the  fluctuation  period  T = 1 //. 

If  we  define  the  apparent  fluctuation  wavelength  as  the  product  of  fluctuation  period 
and  mean  wind  speed,  the  ratio,  F,  of  identification  distance  to  scale  size  can  then  be 
calculated  from  the  mean  slopes  of  Figure  10.  These  ratios  art  ,hown  in  the  figure  and 
range  from  0.24  to  0.50  with  no  apparent  relation  to  height  or  meteorological  para- 
meters. A similar  calculation  was  made  by  Thompson  et  al.  (1968)  who  obtained  a 
ratio  F=  0.4  for  an  array  of  sensors  near  ground  level. 


Fig.  9.  Spatial  correlation  /{( L ) derived  from  maximum  cross  correlations  between  pairs  of  sensors, 
and  the  contribution  from  each  frequency  band  (coherence)  In-line  array  100  m,  IVpq  July  JO  at 

111ft.  with  humidity  sensors. 
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FLUCTUATION  PERIOD,  SECS 

Fig.  10.  Variation  of  identification  distance,  within  which  \'COH>0.6,  with  fluctuation  period. 


FREQUENCY  IN  HERTZ 

Fig.  II.  Phase  of  cross-spectrum  against  frequency  for  various  sensor  spacings.  Data  as  in  Figure  9. 
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Another  way  in  w hich  one  can  examine  the  movement  of  irregularities  is  to  plot  the 
phase  0(f)  of  the  cross-spectrum  against  frequency  for  each  of  the  sensor  pairs; 
since  0 = tan  1 Q/P,  it  is  subject  loan  ambiguity  of  nn  and  this  is  selected  so  as  to  make 
0(f)  continuous.  The  slope  of  the  plot  should  then  lie  related  to  the  apparent  mean 
speed  V of  the  irregularities  as: 


where  d is  the  separation  of  the  appropriate  pair  of  sensors,  f igure  1 1 shows  one  such 
plot  for  four  spacings  and  it  can  be  seen  that  the  points  tend  to  lie  on  straight  lines. 
Substituting  the  mean  slopes  into  Equation  (6),  the  speeds  show  n in  Table  I are  obtained. 
The  average  of  the  four  speeds  is  3.78  m s ' which  agrees  well  with  the  mean  wind 
speed  of  3.5  ms1,  ranging  from  2.5  to  5.0  ms  1 measured  with  a cup  anemometer. 
This  indicates  that  apparent  speeds  derived  from  phase  measurements  give  a useful 
measure  of  the  true  mean  wind  speed. 


TABl.K  1 

Apparent  speeds  derived  from  cross-spectral 
phase  measurements 


Spacing  </  (m) 

1.5 

3.2 

5.75 

9.1 

Speed  V (ms  *) 

3.6 

3.8 

3.8 

3.9 

8.  Conclusions 

We  conclude  by  summarizing  the  main  results  of  this  series  of  measurements,  w hich  in 
some  cases  confirm  and  extend  the  results  of  ground-level  measurements  to  greater 
heights,  and  in  other  cases  provide  new  information  on  atmospheric  parameters  of 
interest  in  the  field  of  radiowave  propagation. 

( 1 ) It  is  possible  to  compute  the  retractive  index  from  fast-response  temperature  and 
humidity  sensors  and  obtain  agreement  with  refractive  index  measurements,  to  10  Hz 
at  least. 

(2)  Generally  the  slopes  of  refractive  index,  humidity  and  temperature  spectra  vary 
over  a wide  range  ( — 1.5  to  - 3.5)  at  frequencies  below  some  I Hz,  and  tend  to  vary  in 
unison.  At  higher  frequencies,  slopes  are  often  close  to  - 5/3  as  predicted  by  the 
theory  of  local  isotropic  turbulence,  and  occasionally  steeper,  up  to  -2.0. 

(3)  It  is  difficult  to  choose  the  scale  over  which  parameters  describing  the  dynamic 
state  of  the  atmosphere  such  as  temperature  gradient,  wind  shear  and  Richardson 
number  should  be  taken  in  order  to  relate  them  to  fluctuations  observed  at  one  point. 
Averages  over  150  m indicate  that  spectrum  slopes  are  close  to  -5/3  when  the 
Richardson  number  is  negative,  the  temperature  gradient  steeper  than-  l.l  C 100  m 1 
and  the  wind  shear  greater  than  1.4  x I0~J  s '. 

Steep  spectrum  slopes  are  frequently  associated  with  stable  atmospheric  conditions 


I 


SMALl  -SCALE  ATMOSPHERIC  STRUCTURE  491 

but  there  is  no  precise  relation  between  slope  and  the  parameters  mentioned  above. 
Thus,  under  stable  conditions  the  parameters  cannot  be  used  to  indicate  the  nature  of 
small-scale  structure. 

There  is  no  apparent  relation  between  spectrum  slope  and  fluctuation  intensity,  but 
intensity  reaches  a maximum  for  a Richardson  number  of  the  order  -0.75. 

(4)  Cross-correlation  measurements  with  vertically-spaced  sensors  show  that  the 
atmosphere  is  sometimes  composed  of  layers  of  order  1-m  in  thickness,  which  on 
occasion  arc  more  marked  in  vapour  pressure  than  in  temperature,  and  sometimes 
give  rise  to  negative  cross-correlations  in  vapour  pressure. 

(5)  Well  above  ground  the  horizontal  spatial  correlation  distances  are  often  large 
and  are  associated  with  steep  spectrum  slopes  at  frequencies  below  some  I Hz. 

(6)  The  cross-correlations  between  sensors  spaced  along  the  mean  wind  direction 
are  larger  than  for  those  set  cross-wind,  horizontally  or  vertically.  Orthogonal  array 
measurements  indicate  large  and  persistent  fluctuations  in  wind  direction  which  may 
be  associated  with  thermals;  this  suggests  taking  caution  in  using  mean  wind  speeds  to 
derive  wave  numbers  from  frequencies  in  spectral  applications. 

(7)  It  is  shown  that,  with  cross-spectral  techniques,  irregularities  can  be  used  as 
tracers  to  obtain  a measurement  of  the  mean  wind  speed  which  agrees  with  direct 
measurements. 

(8)  The  identification  distance,  defined  as  the  distance  within  which  the  coherence 
of  spaced  sensor  measurements  is  greater  than  0.6,  varies  both  in  time  and  space. 
Also  the  ratio  of  identification  distance  to  wavelength,  given  by  the  product  of  wind 
speed  and  fluctuation  period,  varies  over  a range  from  0.25  to  0.5  in  time  and  space 
and  shows  no  simple  relation  to  wind  speed  or  other  meteorological  parameters. 
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